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Modeling statistical properties of motion of a Lagrangian particle advected by a high- 
Reynolds-number flow is of much practical interest and complement traditional studies 
of turbulence made in Eulerian framework. The strong and nonlocal character of La- 
grangian particle coupling due to pressure effects makes the main obstacle to derive tur- 
bulence statistics from the three-dimensional Navier-Stokes equation; motion of a single 
fluid-particle is strongly correlated to that of the other particles. Recent breakthrough 
Lagrangian experiments with high resolution of Kolmogorov scale have motivated grow- 
ing interest to acceleration of a fluid particle. Experimental stationary statistics of La- 
grangian acceleration conditioned on Lagrangian velocity reveals essential dependence 
of the acceleration variance upon the velocity. This is confirmed by direct numerical 
simulations. Lagrangian intermittency is considerably stronger than the Eulerian one. 
Statistics of Lagrangian acceleration depends on Reynolds number. In this review we 
present description of new simple models of Lagrangian acceleration that enable data 
analysis and some advance in phenomenological study of the Lagrangian single-particle 
dynamics. Simple Lagrangian stochastic modeling by Langevin-type dynamical equa- 
tions is one the widely used tools. The models are aimed particularly to describe the 
observed highly non-Gaussian conditional and unconditional acceleration distributions. 
Stochastic one-dimensional toy models capture main features of the observed stationary 
statistics of acceleration. We review various models and focus in a more detail on the 
model which has some deductive support from the Navier-Stokes equation. Comparative 
analysis on the basis of the experimental data and direct numerical simulations is made. 

Keywords: Fully developed turbulence; intermittency; turbulent transport; Lagrangian 
acceleration; conditional acceleration statistics. 

1. Introduction 

In fluid mechanics, acceleration can be defined as the substantive derivative of the 
velocity, 

= "37" = '^tVi + VkdkVi. (1) 
at 

When treated in Eulerian framework, acceleration incorporates the Eulerian local 
acceleration dtVi and the nonlinear advection term Vkduvi which require measure- 
ments of the velocity Vi and temporal and spatial velocity derivatives dtVi and dkVi 
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at fixed point of a flow. Here, dk = d/dx'^ denotes spatial derivative in the labo- 
ratory Cartesian frame of reference, dt — d/dt is time derivative, i. A; = 1, 2, 3, and 
summation over repeated indices is assumed. 

Using Eq. (1) the three-dimensional (3D) Navier-Stokes equation for an incom- 
pressible flow can be written as 

a^ = -p-^d^p + iydlv, + f„ (2) 

where p is the constant fluid density, p is the pressure, v is the kinematic viscosity, Vi 
is the velocity field, and is the forcing, which usually occurs at large characteristic 
spatial scale. 

Measurement of time series Xi (t) of the position of individual tracer particle and 
using a finite-difference scheme allows one to evaluate its Lagrangian velocity Ui{t) 
and acceleration ai{t) as functions of time due to the Lagrangian relations 

Ui{t) = dtx,{t), a,{t) = d'^tx,{t). (3) 

Here, xi = Xi (xofe , t) is the coordinate of an infinitesimal fluid particle viewed as 
a function of the initial position xoi = Xi{0) of the particle and time t G [0, oo). 
With the initial data points xoi (Lagrangian coordinates) running over all the fluid 
particles one gets a Lagrangian description of fluid flow. 

In the Eulerian framework, contributions of the viscous term and the forcing are 
known to be small as compared to that of the pressure gradient term for a certain 
range of scales in a locally isotropic turbulence. Direct analytical evaluation of the 
Lagrangian acceleration ai{xok,t) by using the 3D Navier-Stokes equation for high- 
Reynolds-number (high-Re) turbulent flows is out of reach at present. Thus one is 
led to estimate it theoretically in some fashion. 

Accurate evaluation of the Lagrangian velocity and acceleration in laboratory 
turbulence experiments requires measurement of positions of neutrally buoyant 
tracer particle by using some tracking system to a very high accuracy. One can 
also use direct measurement of the Lagrangian velocity when knowing a precise po- 
sition of tracer particle is not important. In any case, to get Lagrangian acceleration 
one should have experimental access to time scales smaller than the Kolmogorov 
time scale of the flow. One expects that such experimental data give information 
on the pressure gradient term, which is difficult to measure experimentally. 

Growing interest in studying Lagrangian turbulence is motivated by the re- 
cent breakthrough Lagrangian experiment by La Porta, Voth, Crawford. Alexander, 
and Bodenschatz^ (2001), the new data by Crawford, Mordant, Bodenschatz, and 
Reynolds^ (2002), Mordant, Crawford, and Bodenschatz'^ (2003) (optical tracking 
system, R\ = 690, the measured acceleration range is |a|/(a^)^/^ < 60, and is re- 
solved). Mordant, Delour, Leveque, Arneodo, and Pinton^ (2002) (acoustic tracking 
system, Rx = 740, \a\/{a^y^^ < 20, and is not resolved), and direct numerical 
simulations (DNS) of the 3D Navier-Stokes equation by Gotoh, Fukayama, and 
Nakano^ (2002), Kraichnan and Gotoh^ (2003) {Rx = 380, \a\/{aY^^ < 150) and 
Biferale, Boffetta, Celani, Lanotte, and Toschi'^ (2004) {Rx = 280, \a\/{a'^y^'^ < 80, 



stochastic models of Lagrangian acceleration of fluid particle in developed turbulence 3 




Fig. 1. Lin-log plot of the Lagrangian acceleration probability density function P{a) for the 
transverse component of acceleration. Dots: experimental data for the Rx = 690 flow by Craw- 
ford, Mordant, Bodenschatz, and Reynolds.^ Dashed line: the stretched exponential fit (4). The 
acceleration component a is normalized to unit variance. 
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Fig. 2. Lin-lin plot of the central part of the Lagrangian acceleration probability density function 
P{a) for the transverse component of acceleration. Same notation as in Fig. 1. 



is resolved). Experimental results on the 3D Lagrangian acceleration have been 
reported by Mordant, Crawford, and Bodenschatz.^ The classical Reynolds number 
Re is related to the Taylor microscale Reynolds number due to Re — R\/15. These 
experiments and DNS give an important information on single-particle dynamics 
and statistics, and new look to the intcrmittency in high-Rc fluid turbulence. 

The experimental data on a component of Lagrangian acceleration a of 
polystyrene tracer particle in the R\ = 690 water flow generated between two 
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Fig. 3. The contribution to fourth-order moment of Lagrangian acceleration, a^P(a). Same no- 
tation as in Fig. 1. 

counter-rotating disks reveal that the acceleration varies with time in a wild way. 
Statistical description of the acceleration is then used and a huge amount of col- 
lected data have been fitted to a good accuracy by the probability density function 
of a stretched exponential form^'^'^ 

,2 



P{a) = Cexp 



a 



(4) 



il + \bia/b2nbi_ 

Here, hi = 0.513±0.003, 62 = 0.563±0.02, and 63 = 1.600±0.003 are fit parameters, 
and C = 0.733 is the normalization constant. 

Two coordinates z and a;, one along the large-scale symmetry axis and the other 
transverse to it, were measured, while the third coordinate y was taken statisti- 
cally equivalent to the measured transverse component. Lagrangian particles are 
tracked in a small central part of the flow where, in general, high degree of statisti- 
cal isotropy of small scales is expected due to Kolmogorov 1941 (K41) local isotropy 
hypothesis.^ The studied statistically stationary flow is highly anisotropic at large 
scales due to the used specific stirring mechanism. This appears to affect very small 
scales. Namely, one can observe small skewness of the acceleration distribution and 
anisotropy of the acceleration variance, as well as difference in distributions of com- 
ponents of Lagrangian velocity. 

At large acceleration magnitudes, tails of the fitted model distribution (4) decay 
very slowly, asymptotically as exp[— |a|°-^], which implies finitness of the acceleration 
fourth-order moment (a'*) = a'^P{a)da, as confirmed by the experiment.^ The 
flatness factor of the distribution (4) which characterizes widening of its tails (when 
compared to a Gaussian) is F ~ 55.1, which is in agreement with the experimental 
value 

F= {a'^)/{aY = 55±8. (5) 
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Gaussian distribution is characterized by much smaUer value F = 3. 

The Kolmogorov time scale of the R\ = 690 flow is = 0.93 ms. Experimental 
resolution of the scale is very accurate, about 1/65 of r^. Low-pass filtering with 
the 0.23t^ width of the collected 1.7 x 10** data points was used, and the response 
time of optically tracked 46 /xm diameter tracer particle is 0.12r^. Note that certain 
nonzero time scale less than the Kolmogorov time scale is used to derive Lagrangian 
velocity and acceleration values. The width of filter has a limited impact on the 
resulting value of i^, e.g. the use of O.Slr,, width results in about 15% decrease of the 
experimental value of flatness F. Notice that non-ideal response characteristics of 
tracer particle may result in an increase of the effective integral time scale, from the 
Lagrangian integral time scale Tj^ to the Eulcrian integral time scale Te (calculated 
in the comoving frame), as one can show by using Corrsin hypothesis.^" 

The experimental data and stretched exponential fits are shown in Figs. 1, 2 
and 3. One can observe almost symmetric distributions with respect to zero ac- 
celeration and very intermittent character of the Lagrangian acceleration. Namely, 
the pronounced central peak (low accelerations) and long tails (high accelerations) 
make a highly non-Gaussian shape of the acceleration distribution shown in Fig. 1. 
One concludes that the observed fluid-particle dynamics is featured by relatively fre- 
quent acceleration bursts, up to the measured 60 standard deviations. Such extreme 
events occur when the tracer particle is captured by intense small-scale vortical 
structures which are thought to be present in the turbulent flow. These structures 
seem to be distributed randomly in space and time, with large intervals between 
them which are characterized by low-intensity events. As shown by Farge, Pelle- 
grino, and Schneider^^ the most of the turbulent kinetic energy is carried by vortex 
tubes, which are surrounded by a background incoherent flow. 

The long-standing Heisenberg-Yaglom scaling of a component of Lagrangian 
acceleration, ^^'^^ 

(a^) = aou^'/^u-'/^L-^/^, (6) 

was confirmed experimentally^ to a high accuracy, for about seven orders of mag- 
nitude in the acceleration variance, or two orders of the root-mean-square (rms) 
velocity u = [v?)^/"^ (500 < R\ < 980). Here, the Lagrangian velocity u is such that 
the average (u) — 0, oo is the Kolmogorov constant, v is the kinematic viscosity, and 
L is the Lagrangian integral length scale. For R\ < 500 the Heisenberg-Yaglom scal- 
ing was found to be broken. This signals increasing coupling of the acceleration to 
large scales of the flow, and may be related specifically to the large-scale anisotropy 
effect or to "insufficient" developing of the turbulence, or to both of them. 

Very recent experimental data^ on the 3D Lagrangian acceleration in turbulent 
flows with R\ = 285, 485, and 690 show that the three components Ux, ciy, and of 
the acceleration are statistically dependent. For example, the conditional variance 
{a.y\az) increases strongly with the magnitude of ttz- The acceleration magnitude 
\a\ was found to be characterized by the probability density function, which is 
comparable to a log-normal distribution of variance 1 at small and medium values 
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|a|/((|ap) — (|a|)^)^/^ < 25, and by the autocorrelation time of about the integral 
time scale. The autocorrelation time of the direction of the acceleration vector 
a is of the dissipation time scale. The observed two-time-scale character of the 
stochastic dynamics is consistent with previous experimental and DNS results.* 
Assuming the log- normal distribution of the magnitude \a\ and statistical isotropy 
of the acceleration vector a one can straightforwardly derive distribution of each 
component^ 



rametcr. It was shown that the predicted distribution at s = 1 follows experimental 
data points to a good accuracy at small and medium accelerations |aj|/(a?)^/-^ < 25, 
and overestimates them at higher values. The origin of this departure is not clear. 
However, it should be noted that for higher Rx tails of the observed distributions 
become wider, and approach the predicted curve. 

It should be emphasized that the Lagrangian velocity components for the studied 
flow follow Gaussian distribution to a good accuracy. Theoretically, time derivative 
of a dynamical variable does not necessarily follow the same statistical distribu- 
tion as that of the variable. The link between these two sharply distinct distribu- 
tions — Lagrangian velocity and acceleration — can be seen from studying stationary 
statistics of the time increment of a component u of the velocity of individual fluid 
particle. 



For the time scale r of the order of Lagrangian integral time scale Tl (large charac- 
teristic time scale of the flow implied by simple dimensional analysis) the stationary 
distribution of S-ru{t) is approximately of a Gaussian form while for r decreasing 
down to the Kolmogorov time scale (small characteristic time scale of the flow 
implied by simple dimensional analysis) the distribution continuously develops long 
tails, and in a far dissipative subrange it reproduces the acceleration distribution 
shown in Fig. 1. For sufficiently small time scales t < turbulent fluctuations are 
smoothed, and the increment is proportional to the time scale, Sru{t) = Ta(t), to a 
good accuracy. The ratio between the timescalcs, Tl/tjj, is very large for developed 
turbulent flows, and is characterized by Taylor microscale Reynolds number Rx; for 
the studied flow it is of about two orders of magnitude. 

High probability of extreme acceleration magnitudes, as compared to that im- 
plied by the corresponding Gaussian distribution, is associated with the Lagrangian 
turbulence intermittency, which was found to be considerably stronger than the 
Eulerian one. Equivalently, one can say that it is related to an increase of the prob- 
ability to have larger velocity increments in time with decrease of time scale, down 
to the Kolmogorov one (a statistical viewpoint). This is due to the absence of the so 





where erf(x) is the error function, m = i/S/e^**^ for a unit variance and s is a free pa- 



Sru{t) = u{t + r) - u{t). 



(8) 
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called sweeping effect in the Lagrangian frame and the existence of relatively long- 
lived intense vortical structures (vortex tubes) with radii of the order of Kolmogorov 
length T] and total sizes extending up to the integral length scale L. Recent labo- 
ratory experiments by Mouri, Hori, and Kawashima^"' (sec also references therein) 
for boundary layers with Rx = 295-1255 confirm this picture. 

In the traditional Eulerian framework, the isotropic turbulence intermittency is 
understood differently, as an increase of the probability to have larger longitudinal 
velocity differences 

5iv{x) = v{x + /) - v{x) (9) 

on shorter spatial separation scales /, and studied through scaling exponents C^{p) 
of the Eulerian velocity structure functions {{5iv)p) ~ I'' p = 1,2,3,... (a 
structural viewpoint). The velocity difference is taken at the same time instance. 
For Lagrangian velocity structure functions one considers scaling with respect to 
the time scale r, {{Sru)^) ~ r*"^^^. 

Scaling properties of the Eulerian and Lagrangian velocity structure functions 
(statistical moments of 6iv and SrU which characterize their distributions) are tradi- 
tionally used to quantify turbulence intermittency. Here, the intermittency exhibits 
itself as observed deviations from the "normal" scalings predicted by K41 theory. 
The deviations become stronger with increasing order p, with the exception being 
C^(3) = 1 which corresponds to the four-fifth law by Kolmogorov (C(2) = 1 for the 
Lagrangian case).^'^^ 

Under the assumption of balance between the energy injected by driving forces, 
which occur presumably at large spatial scales, and the energy dissipated by viscous 
processes, which are concentrated at small scales, one can restrict consideration by a 
statistically steady state, and focus on intermediate scales, the inertial range, char- 
acteristics of which are universal to some extent for high-Re flows. In the inertial 
range the energy is transferred from large to small scales (the direct energy cascade) 
and viscosity effects are not noticeable. Due to the K41 hypotheses^ certain proper- 
ties of fully developed turbulence in the inertial range are independent on the details 
of initial conditions and forcing (boundary conditions), as well as on the details of 
the energy dissipation. This hypothesis is valid only statistically, in the sense that 
the velocity and acceleration (Heisenberg^^ and Yaglom^^) are viewed as random 
variables. Hence, one is interested in probability density functions and correlators 
of the variables. Complete K41 scale invariancc of the Eulerian velocity difference 
statistics is known to be broken, with the "anomaly" coming from sensitivity of the 
inertial range of scales to large scales of the flow. The breaking occurs also for the 
predicted K41 scalings of the Lagrangian velocity structure functions {{SruY) for 
the same reason. The observed scaling exponents behave as nonlinear functions of 
the order p rather than linear ("normal") ones. This is associated in general with 
the so called dissipative anomaly (mean turbulent kinetic energy dissipation rate 
remains finite when one puts the viscosity parameter to zero) corresponding to a 
strongly nonlinear and non-equilibrium character of a high-Re turbulent system. 
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Possibility to describe high-Re turbulence on the basis of unified statistical princi- 
ples, such as those successfully applied to (quasi-) equilibrium systems, is an open 
problem. The strong and nonlocal character of Lagrangian particle coupling due to 
pressure effects makes the main obstacle to derive turbulence statistics from the 
Navier-Stokes equation; motion of a single fluid-particle is strongly correlated to 
that of the other particles. Studying and accurate modeling Lagrangian dynam- 
ics of a many-particle configuration,^^ and particularly single-particle behavior, is 
currently under elaboration. 

In the inertial range, which extends much for high-Re flows, contribution of the 
pressure gradient to the acceleration variance strongly dominates over that of the 
viscous force. The two contributions do not correlate and the acceleration field can 
be then approximated as the potential one. The DNS pressure gradient data meet 
that of the experimental Lagrangian acceleration.^'^'^ The effects of viscosity and 
intermittency are known to reduce the effective inertial range characterized by the 
ascribed scaling. 

Statistical isotropy and homogeneity of a fully developed turbulent flow at small 
scales are assumed by K41 theory and greatly simplify statistical description of tur- 
bulence. Turbulence intermittency is related to a non-Gaussian statistical behavior 
and is a more subtle matter for theoretical study. Intermittency of the stochastic en- 
ergy dissipation rate at scale I is related to intermittency of dynamics of the system 
that makes a link between the Eulerian and Lagrangian aspects of intermittency. It 
should be emphasized that Eulerian and Lagrangian approaches in studying fluid 
flows are characterized by different theoretical technique and implications, and com- 
pliment each other. 

Modeling statistical turbulence in the Lagrangian frame is important both for 
theoretical implications and applicational studies. Simple models focus on single- 
particle statistical properties and employ Langevin-type equations for one variable. 
Partial justification of the use of one-dimensional models comes from the K41 as- 
sumption on statistical isotropy of the flow at small scales. Away from boundaries 
one can also use the assumption on homogeneity of the flow and discard dependence 
of the Lagrangian tracer on initial position for a sufficiently long evolution. When 
comparing to the experimental or DNS data one usually peaks up one measured 
component of the variable, or uses averaging over all accessible components to get 
higher statistics. Experimentally, a three-dimensional picture is difficult to reach^ 
while DNS naturally gives a full access to it. Some Lagrangian experiments allow 
very accurate resolution of the Kolmogorov time scale with relatively short integra- 
tion time, the others allow to follow individual particle for long time but do not 
resolve the Kolmogorov time scale. DNS is characterized by a high resolution of the 
Kolmogorov time scale and long integration time but at lower Rx achieved because 
of current computational limitations. 

As the experiments and performed DNS (on cubic lattices) do not provide perfect 
isotropy of large-scale forcing (boundary) one naturally expects anisotropy effects at 
all scales. However, these effects are usually small at small scales (even when strong 
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large-scale anisotropy of a high-Re turbulent flow is present) and thus can be ig- 
nored in the first approximation. Influence of anisotropy eff'ects is usually treated 
as a small correction. Local isotropy of a fully developed turbulent flow rather 
strong but very fruitful condition — is one of the main assumptions of K41 phe- 
nomenological theory.^^ However, anisotropy effects at small scales are known to 
be persistent for high Reynolds numbers. For example, persistent anisotropy in the 
skcwncss of velocity derivatives in homogeneous shear flows, which represent one 
of the simplest anisotropic flow, is observed.^* We remind that the forceless 3D 
Navier-Stokes equation is invariant under SO (3) rotational Lie group in coordinate 
space, and symmetry breaking may come only due to the forcing and/or boundary 
conditions. Recently developed SO (3) decomposition theory-'^^ can be used to treat 
anisotropic and isotopic sectors in a rigorous way, and to study how the isotropy 
recovery at small scales happens in Navier-Stokes turbulence. The isotropy recovery 
is partially justified owing to a subleading character of anisotropy found in some 
exactly solvable models. 

Two-time-scale stochastic dynamics in describing the Lagrangian acceleration 
component jointly with the Lagrangian velocity and position of a fluid particle was 
proposed long time ago by Sawford.^° The model equations are 



dta = -{TZ' + T-')a + T£'r-\ + ^2al{Tz' + T^')Tz'T^'L{t), (10) 
dtu = a, dfX = u, (11) 



where 



(12) 



are two time scales, Tl ^ Tr,-, L{t) is zero-mean Gaussian- white noise, Co and ao 
are Lagrangian structure constants, tr^ is the variance of velocity distribution, and 
e is the mean turbulent kinetic energy dissipation rate per unit mass. 

K41 hypotheses of locally isotropic character of high- Re turbulence and similar- 
ity lead to the result that the acceleration field is spatially isotropic, {aiaj) ~ 6ij, 
and the stationary probability distribution of acceleration may depend only on the 
parameters e and v (mean energy flux and viscosity) . The second-order Lagrangian 
velocity structure function {ST-UiS^Uj) ~ Sij also should show spatial isotropy for 
the inertial range of time scales. Thus, a single-component consideration makes a 
sense; u = ui,U2, or us and a = ai,a2, or 03, in laboratory Cartesian frame of 
reference. The constant Co enters the linear scaling of the velocity structure func- 
tion {{dru)"^) = Cqct implied by K41 theory for the inertial range of time scales 
Tr, <^ T <^ Tl- Since the form of this two-time correlator is similar to that of 
the displacement of usual Brownian particle, the velocity of fluid particle in the 
inertial range can be thought of as a Brownian-like motion with the "diffusion" 
coefficient Coe. In other words, the velocity is a stationary stochastic process with 
independent increments. For time scales smaller than the Kolmogorov time scale 
Tr, the predicted scaling is quite different, (((J^u)^) = aoe'^^^i'~^^^T'^, and directly 
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corresponds to Heisenberg-Yaglom scaling law (6); e = u'^/L, rj = {v^/iy^^, and 
R\ = {Ib/ivY^'^u'^. This theory also predicts e/r decay of the autocorrelation of 
acceleration component {a{t)a{t + t)) , when one imposes its independence on i/, for 
time scales t much bigger than r,;. For the velocity autocorrelation {u{t)u{t + r)) 
the prediction is that it decays considerably only for r of the order of Lagrangian 
integral time scale Tj,. The uncorrelated character of Lagrangian acceleration then 
could be used to build a first approximation for time scales within the inertial range. 
Note however that the predicted decay of {a{t)a(t + r)) across this range is rather 
slow. Qualitatively K41 relations mean that the components of Lagrangian acceler- 
ation and velocity are associated mainly with small and large scales of a developed 
turbulent flow respectively. 

Sawford model (10) predicts stationary Gaussian distributions for both the ac- 
celeration and velocity reflecting the used uncorrelated character of fluctuations 
and is consistent with K41 picture. One of the extensions of this model is due to 
replacement of e by stochastic energy dissipation rate e, and assuming that it is 
lognormally distributed in correspondence to the refined Kolmogorov 1962 (K62) 
approach. Such extensions can be used to fit the observed highly non-Gaussian 
shape of the acceleration distribution shown in Fig. 1. 

Recent Lagrangian experiments and DNS by Mordant, Delour, Leveque, Ar- 
neodo, and Pinton* and DNS by Chevillard, Roux, Leveque, Mordant, Pinton, and 
Arneodo^^ show that certain long-time correlations and the occurrence of very large 
fluctuations at small scales dominate the motion of a fluid particle, and this leads 
to a new dynamical picture of turbulence. This requires effective models on how to 
account for the specific long-time correlations along a particle trajectory which are 
viewed as a key to intermittency in turbulence. 

While it is evident that the 3D Navier-Stokes equation with a Gaussian-white 
random forcing belongs to a class of non-linear stochastic dynamical equations for 
the velocity field with which one can associate some generalized Fokker-Planck 
equations or apply a path-integral method, it is a theoretical challenge to make 
a link between the Navier-Stokes equation (2) and phenomenological stochastic 
models of Lagrangian acceleration. 

Recent approach by Friedrich^'' can be traced back to the so called probability 
density function method originated by Oboukhov^^ and developed by Pope^® and 
Sawford'^^ (see also references therein). Friedrich has shown that one can obtain 
infinite chain of evolution equations for joint Lagrangian n-point probability den- 
sity functions and closed equation for the associated probability density functional 
which stem from the incompressible 3D Navier-Stokes equation in the Lagrangian 
frame; see also work by Heppe.^* Evolution equation for the single-particle dis- 
tribution function f{v,r,Xo,t) = {6{v — u{xo,t))6{r — x{xo,t))), where u{xo,t) 
and x{xo,t) are Lagrangian velocity and position respectively, includes integral of 
pressure gradient and dissipation operators acting on mixed Eulerian-Lagrangian 
equal-time two-particle distribution function, and so on. Particularly, he derived a 
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generalized Fokker-Planck equation (with memory term) for a single-particle proba- 
bility distribution of Lagrangian velocity increments by using certain closure scheme 
partially justified for high-Rc homogeneous isotropic turbulence. The approach nat- 
urally leads to consideration of acceleration covariances conditional on Lagrangian 
velocity and position which correspond to a three-point distribution function. Such 
a conditional dependence was dropped in order to reduce the Fokker-Planck equa- 
tion, which nevertheless accoimts for time intc;grated efFcK;ts. This approximation 
means particularly that the correlation between acceleration fields at space-time 
points {r,t) and (r — l,t') does not depend on the velocity of a fluid particle at 
(r — I, t'), where I = v{t — t'). K41 theory is used to derive general form of the two- 
point two-time acceleration autocorrelation function, which approximates diffusion 
term, for the inertial range, whereas the drift term vanishes identically because of 
the ignored conditional dependence. Power-law form for unknown function entering 
the diffusion term in the Fokker-Planck equation for modulus of velocity was used, 
with the exponent being treated as a free parameter. This leads to consideration of 
a class of continuous-time random walk of the velocity featured by non-Markovian 
behavior, which is contrasted to Markovian treatment (no memory effects) under- 
lying Oboukhov model^"'' with Gaussian distributed Lagrangian acceleration. The 
resulting equation is analytically tractable, and its solution is presented in the form 
of definite integral. Timescale dependence of the free parameter was used to fit the 
experimental data on statistics of Lagrangian velocity increments in a wide range of 
timescalcs.^^ The introduced timescale dependence requires a justification since this 
parameter was treated constant when solving the evolution equation. The closure 
scheme provides the following scaling behavior of the Lagrangian velocity distribu- 
tion: P{u, t) = t~'^^'^Ps{ut~^/^). Importance of this approach is that it has deductive 
support from the Navier-Stokes turbulence and directly accounts for the memory 
effects. 

Lagrangian acceleration statistics conditional on the same component of 
Lagrangian velocity was studied experimentally by Mordant, Crawford, and 
Bodenschatz.^ These data add a very useful information on the Lagrangian in- 
termittency as well as allow one to check implications of refined stochastic models, 
which describe distribution of the acceleration conditional on velocity. 

The conditional acceleration probability density function P{a\u) at a set of fixed 
velocities u ranging from to 3.1 (in rms units) was found to be of approximately 
the same stretched exponential shape as that of the unconditional acceleration P{a) 
shown in Fig. 1. Theoretically, the distribution P{a) can be calculated with the use 
of P{a\u) by integrating out u in P{a\u) with some (independent) distribution of 
u. The experimental conditional acceleration variance {a^lu) was found to increase 
in a nonlinear way with the increase of magnitude of velocity u. Dependence of the 
acceleration variance on the velocity magnitude breaks local homogeneity of the 
flow assumed by K41 theory, and is a prerequisite to describe turbulence intermit- 
tency. One therefore should admit influence of larger scales when describing the 
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small-scale dynamics by supposing that the intense structures are characterized by 
both the large and small time scales in the Lagrangian framework. The conditional 
mean acceleration {a\u) was found to be nonzero and increases for higher veloc- 
ity magnitudes that reflects the large-scale anisotropy effect of the studied flow. 
Recent DNS result by Biferale, Boffetta, Celani, Devenish, Lanotte, and Toschi^^ 
based on the analysis of 3.6 x 10® data points also shows an essential dependence 
of the acceleration variance on magnitude of large velocities. These findings are 
consistent with the understanding that the long-time correlations along a particle 
trajectory dominate the motion since Lagrangian velocity is characterized by the 
"energy-containing" scales of a turbulent flow. 

The aim of the present paper is to review simple Langcvin-type single-particle 
modeling approach and make a comparative analysis of different recent models of 
Lagrangian acceleration, on the basis of recent Lagrangian experimental data and 
direct numerical simulations of high-Rc isotropic turbulence. We restrict consid- 
eration by steady-state Lagrangian single-particle statistics. Most of the reviewed 
models are one-dimensional. Such models can shed some light to a more realistic 
three-dimensional modeling. We also briefly review some recent results of alterna- 
tive approaches, — multifractal description and multifractal random walk model of 
homogeneous and isotropic turbulence — to provide the reader with a current view 
on the problem. 

The layout of the paper is as follows. 

In Sec. 2 we briefly review recent multifractal approaches to the Lagrangian and 
Eulerian intermittency. The formalism by Chevillard, Roux, Leveque, Mordant, Pin- 
ton, and Arneodo^^ is a Lagrangian version of Eulerian multifractal approach and 

describes statistics of Lagrangian velocity increments in a wide range of time scales, 
from the integral to dissipative one. Fine structure of the range of spatial scales 
smaller than the inertial range has been considered by Chevillard, Castaing, and 
Leveque. ^° Arimitsu and Arimitsu"^^ have constructed multifractal cascade model 
to derive Lagrangian acceleration distribution by making a link between cascade 
picture of isotropic turbulence and Tsallis nonextensive statistics formalism. "^^ 

In Sec. 3 we review some recent one-dimensional Langevin-type models of the 
Lagrangian acceleration in developed turbulence. In Sec. 3.1, we outline implica- 
tions of the models by Beck'^^''^^ with the underlying ^-square (Sec. 3.1.1) and 
log-normal (Sec. 3.1.2) distributions of the model parameter /3, and the x-square 
Gaussian model. We review results of the so called Random Intensity of Noise 
(RIN) approach'^^ to specify the probability density function /(/3) which is based 
upon relating /3 to normally distributed velocity u (Sec. 3.1.3). This formalism en- 
ables to reproduce x-square and log-normal distributions of /? as particular cases. 

A nonlinear Langevin and the associated Fokkcr-Planck equations obtained 
by a direct requirement that the probability distribution satisfies some model- 
independent scaling relation have been recently proposed by Hnat, Chapman, and 
Rowlands^^ to describe the measured time series of the solar wind bulk plasma 
parameters. We find this result relevant to fluid turbulence since it is based on a 
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stochastic dynamical framework and leads to the stationary distribution with expo- 
nentially truncated power law tails, similar to that obtained in the above mentioned 
RIN models. This model is reviewed in Sec. 3.2. 

Recent second-order and third-order Langevin stochastic models of Lagrangian 
acceleration developed by Reynolds'^®''^^''"'''*^ are reviewed in Sec. 3.3. The second- 
order model generalizes Sawford stochastic model (10) while the third-order model 
introduces hypcr-acceleration (substantive derivative of the acceleration) and the 
associated time scale. When neglecting third-order processes one recovers a second- 
order model. Reynolds-number effects are incorporated into the second-order model, 
which is applicable at large time scale. Such a modeling of accelerations in homoge- 
neous anisotropic turbulence has been recently made by Reynolds, Yeo and Lce."*^ 
Reynolds and Veneziani'*'^ have shown importance of trajectory-rotations and that 
non-zero mean rotations are associated with suppressed rates of turbulent dispersion 
and oscillatory Lagrangian velocity autocorrelation functions. Particularly, due to 
the developed extended second-order model, non-zero conditional mean acceleration 
endows trajectories with a preferred sense of rotation. 

The Navier-Stokes equation based approach to describe statistical properties 
of small-scale velocity increments, both in the Eulerian and Lagrangian frames, 
was developed in much detail by Laval, DubruUe, and Nazarenko;^^ see also recent 
work by Laval, Dubrulle, and McWilliams.^^ This approach introduces nonlocal 
interactions between well separated scales, the so called elongated triads, and is 
referred to as the Rapid Distortion Theory (RDT) approach. This approach is con- 
trasted with Gledzer-Ohkitani-Yamada (GOY) shell model, in which interactions 
of a shell of wave numbers with only its nearest and next-nearest shells are taken 
into accoimt. In Sec. 3.4 we outline results of this approach and focus on the pro- 
posed one-dimensional Langevin-type model of Lagrangian small-scale turbulence 
to which we refer as Laval-Dubrulle-Nazarenko (LDN) model. This model includes 
Gaussian-white additive and multiplicative noises with constant intensities, while 
local interactions are accounted for by introducing a turbulent viscosity. LDN-type 
model for Lagrangian acceleration exploits such a simple form of the noises, which 
represent effects of stochastic distortion produced by large scales. 

In Sec. 4 we represent in some detail qualitative and quantitative comparative 
analysis of the one-dimensional LDN-type model at zero correlation between the 
noises and simple RIN models.*^ 

In Sec. 5 we review very recent models of the conditional acceleration statis- 
tics by Sawford, Ycung, Borgas, Vedula, La Porta, Crawford, and Bodcnschatz,"*^ 
Reynolds,^" and Biferale et al.'^'^ We present our study''^ on the conditional prob- 
ability density function P{a\u) where the one-dimensional LDN-type model with 
mutually (5-correlated Gaussian-white additive and multiplicative noises is taken as 
a constitutive model and certain model parameters are assumed to depend on the 
amplitude of Lagrangian velocity u. We also present results of a complete quantita- 
tive description of the available experimental data^'^'^ on conditional and uncondi- 
tional acceleration statistics within the framework of a single LDN-type model with 
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a single set of fit parameters. 

In Sec. 6 we briefly review recent results of the application of multifractal random 
walk theory by Muzy and Bacry^^'^*^ to developed turbulence. This approach allows 

one to go beyond modeling of the Lagrangian velocity of fluid particle by Gaussian 
process to include Poisson process, and to use Kolmogorov-Levy-Khinchin theory 
of stochastic processes with independent increments. 

2. Multifractal approaches 

Recently, Chevillard, Roux, Lcvoquc, Mordant, Pinton, and Arneodo^^ have con- 
structed an appropriately rccastcd multifractal approach, which is widely used in 
Eulerian studies of turbulence, to describe statistics of Lagrangian velocity incre- 
ments in a wide range of timescales, from the integral to dissipative one. The re- 
sulting theoretical distribution reproduces continuous widening of the velocity in- 
crement probability density function (PDF) with the decrease of time scale, from 
a Gaussian-shaped to the stretched exponential, as observed in Lagrangian exper- 
iments carried out at Gornell-^'^'^ and ENS-Lyon,^'^^ and DNS of the 3D Navicr- 
Stokes equation. Two global parameters (Reynolds number and Lagrangian integral 
time scale) and two local parameters (smoothing parameter and intermittency pa- 
rameter) with a parabolic singularity spectrum were used to cover the data in the 
entire range of time scales. 

At dissipative time scale the obtained PDF flts the experimental data on La- 
grangian acceleration to a good accuracy. The cumulant analysis made in this ap- 
proach provides an understanding of the observed departures from the scaling when 
going from the integral to dissipative time scale. The used parabolic singularity 
spectrum D{h) is a hallmark of the log-normal (Kolmogorov 1962) statistics and 
reproduces well the left-hand side (corresponding to intense velocity increments) of 
the observed curve, which is centered at 0.58 (> 1/2), but increasingly deviates at 
the right-hand side (rhs) of it (corresponding to weak velocity increments). Another 
widely used statistics, the log-Poisson one, was shown to imply departure from the 
Lagrangian observations in the same manner. The acceleration statistics conditional 
on velocity was not considered in this work. 

The basic assumption of the multifractal approach is to relate Lagrangian ve- 
locity increments at different time scales to each other, 

u{t + t)- u{t) = I3{T/T){u{t + T)- u(t)), (13) 

by using independent random function /3{t/T), where the time scale r is such that 
T < T and T is fixed at the order of Lagrangian integral time scale T^. This 
relation is understood as a statistical law. When considering the function (3{t /T) 
to be deterministic, one ends up with a monofractal (monoscale, or self-similar) 
picture, well-known example of which is Brownian motion. Random character of 
/3(r/T) leads to a multiscale behavior^^'^^ of the stochastic velocity, for which scaling 
properties of structure functions can be readily derived. ^'^ The scaling exponent C(p) 
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of the p-order structure function is linear for monofractal processes and nonlinear 
in p for multifractal processes. The simplest example of multifractal process is given 
by log-normal distribution of (3, for which case the resulting scaling exponent is a 
second-order polynomial in p (parabola). 

The multifractal approach has been managed to describe both the inertial and 
dissipative range of time scales. Thus the effect of dissipation has been accounted. 
The condition that the local Reynolds number at some time scale between the in- 
ertial and dissipative ranges is of the order of unity was used, and the Lagrangian 
integral time scale and Reynolds number, which are available from the experimen- 
tal data, are explicitly incorporated into the formalism. The remaining two free 
parameters were used for fitting. The obtained PDF of the Lagrangian velocity 
increments is symmetric and includes both the distinct regimes in a unified way 
by using Batchelor's interpolation formula, which contains the smoothing parame- 
ter. The other free parameter is the Lagrangian intcrmittcncy parameter measuring 
second-order nonlinearity of the scaling exponent of the Lagrangian velocity struc- 
ture function. Gaussian shape for the PDF of velocity increments at the integral 
time scale was used, and the Kolmogorov scaling of the second-order Lagrangian 
velocity structure function {{Sru)"^) ~ r was assumed for the inc^rtial range. The 
results of this approach are in a good agreement with various sets of experimental 
and DNS data. 

In a very recent paper Chcvillard, Castaing, and Leveque^" have considered a fine 
structure of the range of scales smaller than the inertial range, within the framework 
of Eulerian multifractal approach. Below, we briefly review and discuss results of 
this work as it concerns very small scales with which fluid-particle acceleration 
is generally associated. Indeed, Hcisenbcrg-Yaglom scaling law (6) shows that the 
acceleration essentially depends on the viscosity v so that it is associated with very 
small scales Z < of a turbulent flow for which viscous effects are known to dominate 
over inertial effects. 

They introduced the so called far-dissipative range of spatial scales, I < -q^, and 
the near-dissipative range, I G [77_,?7_|-], and found that the Eulerian intermittency 
grows faster across the scales in the near-dissipative range as compared with that 
in the inertial range, Z > 77+, with decreasing separation length I of the longitudinal 
velocity difference; the Kolmogorov length scale 77 is such that r]- < rj < 1]+. This 
observed phenomenon has been qualitatively described and attributed to the rein- 
forcement of contrast between intense and quiescent motions due to the gradually 
increasing scale-dependent viscous cutoff effect when going to smaller scales, from 
ry+ to rj- . As the typical scalings in the inertial and far-dissipative ranges are known 
one can compute relationship between the values of cumulants at the endpoints of 
the near-dissipative range. The so called 9/4 amplification law for the intcrmittcncy 
has been established: the parameter measuring intermittency increases more rapidly 
during the crossover from the inertial range to the far-dissipative range as compared 
to that in the inertial range. Highly remarkably, this result has been found to be 
independent on the Reynolds number. 
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In general, the far-dissipative range is characterized by a strong viscous damping 
and eventual saturation of the intermittency with decreasing scale which reaches 
some highest possible value for a given turbulent flow. Reynolds number depen- 
dence of the left and right endpoints rj^ and 77+ has been established and described: 
the near-dissipative range of scales becomes wider for higher Reynolds numbers 
exhibiting approximately VlnRe dependence but the inertial range grows faster, 
approximately as In Re. Scaling properties of Eulcrian velocity structure functions 
were found to be non-universal in the near-dissipative range of scales as they depend 
on the Reynolds number and the order of structure function. The experiment and 
DNS were made for Re = 11750 and Re = 1070 flows respectively'"^ for which pro- 
nounced near-dissipation range is observed. For sufficiently high Reynolds numbers 
the near-dissipative range can be taken negligible when compared to the inertial 
range, and the Kolmogorov length 77 ~ ?7_ ~ 77+ becomes the only scale character- 
izing dissipativc range that appears to bo in accord to the original K41 theory. 

While it is evident that viscous effects are ultimately responsible for such a 
behavior of the intermittency parameter in the near-dissipative range and that the 
Kolmogorov scale can be defined by setting local Reynolds number to unity, it is 
a rather subtle matter to relate corrections to intermittency coming from viscous 
cutoff scales with Navier-Stokes turbulence dynamics. Below we present a tentative 
picture. 

In general, one naturally expects a specific crossover from the mechanism of 
downscale energy-transfer operating within the inertial range (the direct energy 
cascade with negligible effect of viscosity) to the mechanism of strong viscous dissi- 
pation of energy at very small scales of the flow. These two mechanisms are evidently 
of a quite different character, and matching one to the other requires specific inter- 
mittent structures which span the near-dissipative range. 

Possible downscale phenomenological picture is that small scales become much 
more intense as one riches the Kolmogorov length scale; r/. whic;h is a characteristic 
size of very intense vortical structures responsible for intermittent bursts in this 
range. The viscous damping tends to weaken and destroy such intense vortices and 
their correlations, and terminate formation of smaller-radius vortex tubes for a 
given Reynolds number. For local Reynolds numbers less than unity one therefore 
expects no cascade picture except for a single (the smallest-radius) vortex tube, 
while for that bigger than unity cluster of highly correlated vortex tubes is likely 
to be present. The 9/4 amplification law might be due to relative under-population 
and/or more coherent character of such intense vortical structures as compared to 
that of vortical structures in the inertial range. 

Vortex tubes are viewed as the most elementary structures in turbulence. 
Recent rough-wall boundary-layer experimental analysis by Mouri. Hori, and 
Kawashima^'* of spatial distribution of the small-scale vortex tubes with character- 
istic radii of 6.1r]-7Ari shows that the probability to find small separations between 
the tubes is considerably higher than that of large separations. The experimental 
data are due to one-dimensional cut. The vortex tubes have been identified using 
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enhancements of the velocity increment above certain threshold, which eliminates 
detection of weak tubes and other low-intensity structures. The data give access 
down to fractions of Taylor microscale A = [2{v'^) / {{dxv)'^)]^/'^ of the studied wind 
tunnel flows. The length parameter A increases from 3477 to TOry for the flows with 
Rx from 295 to 1255, while Kolmogorov scale rj decreases from 0.06 cm to 0.01 cm 
respectively; R\ = {v'^y^'^X/v. The radius of vortex tube is identified by the po- 
sition of maximum of the obtained Burgcrs-likc antisymmetrizcd velocity profile 
with respect to zero. The finding is that large separations are distributed in a ran- 
dom and independent way due to the observed exponential tail of the distribution 
whereas smaller separations, below few Taylor lengths, occur increasingly frequent. 
Superposition of two exponential functions was used as a model which fits the ex- 
perimental distribution at large and middle separations (from 25A to 5A) to a good 
accuracy but increasingly underestimates it at smaller separations (less than 5A). 
This indicates that the vortex tubes tend to cluster together and correlate to each 
other at the middle and small separation scales. Shapes of the experimental spa- 
tial distributions, argument of which is expressed in Taylor microscale units, are 
found approximately the same for a wide range of Reynolds numbers. This can be 
treated as a universal character of the spatial distribution of vortex tubes of the 
same small dimensionless radius (in Kolmogorov length units) in high- Re turbulent 
flows, Rx > 400. Reynolds-number dependence of the tube parameters was studied 
as well: the radius scales linearly with Kolmogorov length, i?o ^ 77, and the max- 
imum circulation velocity scales linearly with the rms velocity, Vq ~ (f^)^/^, for 
400 < Rx < 1255; the linear scaling of the velocity with respect to the Kolmogorov 
velocity Uk = (i^e)^/^ which is acceptable within the framework of K41 theory is 
not observed. 

As the tendency is that small-radius vortex tubes form clusters and correlate to 
each other at small separations (which incorporate the inertial range of scales), one 
expects more sparse character of vortical structures at smaller scales. Existence of 
vortex tubes with smaller radii is supported by the observation that the Reynolds 

1 /2 

number Reg — VqRq/v characterizing circulations of the vortex tubes scales as R^ 
(i.e. not constant), with the value increasing from Reo = 32 to 62 for Rx = 295 to 
1255. Despite that high Rcq may result in less stable character of the vortex tubes 
they have a rather long lifetime, which is of the order of large-eddy turnover time 
L/(t'^)^/^. This effect may be due to the clustering. 

Since one of the characteristic parameters of the small-scale vortical structures 
essentially depends on Rx it is then natural to observe Reynolds-number dependence 
of the Lagrangian acceleration statistics;^ the level of intermittency increases with 
Reynolds number. 

We note that for vortex tubes with small Rcq < 1 one expects no clustering. High 
probability to find the same-radius vortex tubes with Rcq > 1 separated by small 
distance suggests that they form coupled pairs, or multipole clusters consisting 
of several vortex elements in a more general case. Conservation of the enstrophy 
moment for a vortex pair can be used to find minimal distance between centers 
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of vortex tubes which is estimated to be of about 3.1-3.7 radius of the tube. 
Hence the probabiUty of separations smaller than about 3 radius 21r] for the 
studied tubes) should tend to zero. The experimental distribution^^ is currently 
available for separations down to few radii for which a tendency to reach maximum 
and drop with decreasing separation is indeed observed. Accurate resolution would 
allow to identify the relative separation corresponding to maximum probability and 
study its Reynolds-number dependence. This could help to investigate the vortex 
clusters which appear to be typical small-scale objects in high-Re turbulent flows. 
The vortical structures are advected by a noisy background flow, such as the eddy- 
noise forcing^^ and the energy-equipartition based random forcing.^^ 

It is important to note that just higher mean intensity of small scales does 
not guarantee increase of intermittency since the latter is associated mainly with 
(nonlocal) interactions of the small scales with much larger scales, as argued recently 
by Laval, DubruUe, and Nazarenko.**^ 

The specific increase of intermittency in the near-dissipative range indicates 
just a noticeable character rather than the overall essential gradual increase of the 
effective viscous damping across the near-dissipative range. This is contrasted to 
what happens in the far-dissipative range which corresponds mainly to the interior 
and surrounding of small-scale intense vortex tubes where strong viscous dissipa- 
tion occurs. However, one should keep in mind that some other parameters control 
intermittency as well when one treats turbulence intermittency as a nonlocal phe- 
nomenon. For example, (a) higher intensity of stochastic large-scale strain coupled 
to small scales and (b) stronger large-scale correlations across the inertial range 
produce considerable increase of the intermittency at the small scales. 

It seems that more sparse character of coherent structures at small scales, rather 
than some more or less intense local interactions, could be directly associated with 
the rapid increase of intermittency in the near-dissipative range. The role of viscous 
damping is essential here, but it is likely restricted to the strong damping effect at 
the smaller-scale end of the near-dissipative range. In general, this could be viewed 
as a phenomenon related to maintaining of the mean energy flux downscale to the 
far-dissipative range. It seems that the local Reynolds number varies much due 
to high inhomogeneity of the flow at the small scales so that more refined tools as 
compared to the usual Fourier transform are required to capture details of the small 
scales. 

More detailed analysis and interpretation of experiments resolving Kolmogorov 

scale are required for the near-dissipative range which is currcuitly much less under- 
stood than the inertial range. In particular, whether the flatness factor of the distri- 
bution of Lagrangian velocity increments exhibits a pronounced range of timescales 
similar to the near-dissipative range of spatial scales is still an open question. "^'^^ 

In a very recent paper Bifcrale, Boffetta, Celani, Lanotte, and Toschi*" have 
presented interesting results of DNS of Lagrangian transport in homogeneous and 
isotropic turbulence with Rx up to 280, a very accurate resolution of dissipative 
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scale, and an integration time of about Lagrangian integral time scale. This is con- 
trasted to the experimental optical tracking^'^ which gives access to the resolution 
of about 1/65 of the Kolmogorov time scale and integration time of about few 
Kolmogorov time scales of the R\ = 690 flow, and to the experimental acoustic 
tracking'* which enables the integration time of the order of Lagrangian integral 
time scale but no access to time scales less than the Kolmogorov time scale of the 
R\ = 740 flow. While the value Rx = 280 is not as high as those in the exper- 
iments, an additional advantage of the DNS data is that it gives access to a full 
three-dimensional picture of the flow and high statistics is reached by "seeding" and 
tracking simultaneously millions of Lagrangian tracers. 

In the subsequent work, Bifcralc, Boffctta, Cclani, Dcvcnish, Lanotte, and 
Toschi^^ have shown how the multifractal formalism offers an alternative approach 
which is rooted in the phenomenology of turbulence. The Lagrangian statistics was 
derived from the Eulcrian statistics without introducing ad hoc hypotheses. Shc- 
Leveque empirical formula for the Eulerian scaling exponents has been used and 
time scale is related to the length scale by using the assumption that Eulerian 
velocity diflEerences are proportional to Lagrangian velocity increments, Siv ~ ^r^^- 
Although the formalism is not capable to account for small acceleration values (typi- 
cal situation for the multifractal approach), the obtained acceleration PDF captures 
the DNS data to a good accuracy in the tails, with acceleration values ranging from 
about \a\/ (o,^)^/^ = 1 up to 80. Alas, one can observe an overestimation in this range 
which can be clearly seen from the predicted contribution to fourth-order moment, 
a'^P{a), as compared to the DNS data. High degree of isotropy of the simulated 
stationary flow suggests statistical equivalence of Cartesian components of acceler- 
ation aligned to fixed directions, and the resulting DNS acceleration distribution 
obtained by averaging over the components has been found, as expected, with no 
observable asymmetry with respect to a — > —a. 

Also, acceleration variance conditional on the velocity has been derived^^ within 
the same multifractal approach and compared to the DNS data. We will consider 
this issue below in Sec. 5. 

Recent multifractal cascade model by Arimitsu and Arimitsu^^ implies La- 
grangian acceleration distribution, which fits DNS acceleration data'''"^ to a very 
good accuracy. The model is based on the analysis of scale invariance of the 3D 
Navier-Stokes equation for high Reynolds numbers, and on the assumption that 
singularities due to the invariance distribute themselves multifractally in physical 
space. The guiding principle is an extremum of Tsallis nonextensive entropy^^ under 
certain constraints from which distribution function P{a) of singularity exponent a 
is obtained. Basically, two fit parameters determine the resulting distribution: the 
total number of "steps" in turbulent cascade, n, and the intcrmittcncy exponent fi. 
The ascribed "eddy size" decrement factor is 2. Each step assumes statistical inde- 
pendence of the corresponding flow modes within the multifractal scaling range. The 
acceleration statistics is obtained from the scaling Spm/Spo = [Im/loY"^^ , where 
5pm is the pressure difference across the separation length 1^ = 2~'"Zo and Iq is the 
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turbulence integral scale. At the step m = n the cascade is terminated and one ex- 
pects good approximation for the pressure gradient. Minimum and maximum values 
of a for the singularity spectrum /(a) are related to Tsallis entropic index q by 

(1 - = l/a_ - l/a+, (14) 

f{ct±) = 0. The representation for spectrum corresponding to the cascade model of 
isotropic turbulence is found as follows: 

/(a) = 1 + (1 - log2[l - (a - ao)y{Aa)% (15) 

where (Aa)'^ = 2X/[{1 - (?)ln2] and q, X, and ao arc determined from /i; a± = 
ao ± Aa. Energy conservation and definition of the intermittency exponent were 
used. 

It was argued that the acceleration distribution should include two parts: one 
coming from the above multifractal analysis and the other corresponding to contri- 
bution of dissipative term (the so called "thermal fluctuations" and/or measurement 
errors). The first part determines shape of the tails whereas the second part makes 

the core of distribution (small acceleration magnitudes) with another parameter q' 
entering model Tsallis distribution. Very good fit to the DNS acceleration data is 
obtained for the values n = 18, /x = 0.240 (q = 0.391), and q' = 1.7. However, as 
pointed out by Kraichnan and Gotoh,^ the total number of steps for the simulated 
R\ = 380 flow should not exceed n = 9 to provide consistent treatment of the 
cascade, and that there is no way to flt the tails at n = 9 by tuning the value of jj,. 

3. One-dimensional Langevin toy models of Lagrangian 
acceleration in turbulence 

In this Section, we outline results of some recent one-dimensional Langcvin-type 
models of Lagrangian fluid particle acceleration in a developed turbulent flow. 

Modeling of the Lagrangian acceleration dynamics can be naturally made by 
employing Langevin-type equation, which contains time derivative of the accelera- 
tion, so that random forces responsible for the time evolution of acceleration of a 
fluid particle are related to the time derivative of the rhs of Eq. (2) treated in the 
Lagrangian frame. 

Various one-dimensional models were suggested recently to describe Lagrangian 
acceleration statistics among which we mention the x-square modeP'''^^ and 
the log-normal model"^** by Beck which are based on the Tsallis nonextensive 
statistics^^ inspired approach,^^'^^'^*^'^^ the second-order and third-order models 
by Reynolds'*^ '''^'''"'^^ which extends the model by Sawford,^" the x-square Gaus- 
sian model, ■^^'^^ and the model with underlying normally distributed Lagrangian 
velocity fluctuations.^^ 

It is worthwhile to note that the idea to use stochastic averaging over random 
variance of intermittent variable was used long time ago by Castaing, Gagne, and 
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Hopfinger.^^ Their propagator approach is related to the so called Markovian de- 
scription by Renner, Peinke, and Priedrich^^ as shown by Donkov, Donkov, and 
Grancharova;^^ see also work by Amblard and Drossier. 

Review and critical analysis of the applications of various recent noncxtcnsive 
statistics based models to the turbulence have been made by Gotoh and Kraichnan.^ 
An emphasis was made that some models lack justification of a fit from turbulence 
dynamics although being able to reproduce experimental and DNS data to more or 
less accuracy. Deductive support from the 3D Navier-Stokes equation was stressed to 
be essential for the fitting procedure to be considered meaningful; see also Ref. 46 for 
a review. Also, Zanette and Montemurro^^ have argued recently that the connection 
between specific non-thcrmodynamical processes and non-extensive mechanisms is 
generally not well defined. 

3.1. Simple RIN models 

Tsallis nonextensive statistics^^ inspired formalism^^'^^'^^ was recently used by 

Beck^^'^^ to describe Lagrangian statistical properties of developed turbulence; see 
also Rcfs. 38, 57, 59. In recent papers, ^^'^^'^^ we have made some refinements of 
this approach. 

In this approach, PDF of a component of acceleration of infinitesimal fluid par- 
ticle in developed turbulent flow is found due to the equation 

/■oo 

P{a) = / dp P{a\/3)m, (16) 
Jo 

where P{a\P) is PDF of a conditional on fi. This distribution is associated with 
a surrogate dynamical equation, the one- dimensional Langevin equation for the 
Lagrangian acceleration a, 

dta = ^F{a)+aL{t). (17) 

Here, dt denotes time derivative, F{a) is the deterministic drift force, 7 is the drift 

coefficient, measures intensity of the noise, a strength of the additive stochastic 
force, and L{t) is Gaussian-white noise with zero mean, 

{L{t)L{t')) = 25{t-t'), {L{t))=0. (18) 

The averaging is made over ensemble realizations. Short-time correlated force, which 
is approximated here by L{t), is assumed to come as a combined dynamical effect of 
the flow modes. This force can be viewed as a "background" stochastic force which 
acts along a particle trajectory. 

For constant model parameters 7 and a, this usual Langevin model ensures that 
the stochastic process a{t) defined by Eq. (17) is Markovian. The PDF P{a\(3) of 
the acceleration at fixed /3 = 7/0-^ can be found as a stationary solution of the 
corresponding Fokker-Planck equation 

dtP{a, t) = a„[-7F(a) + cj''d^]P{a, t), (19) 
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where da = d/da. This equation can be derived from the Langevin equation (17) 
using the noise (18) either in Stratonovich or Ito interpretations. Particularly, for a 
linear drift force F{a) = —a, the stationary PDF, i.e. dtP{a, t) = 0, is of a Gaussian 
form, 

P(a|/3) = C(/3)exp[-/3aV2], (20) 

where C(/3) = (3/ (27r) is a normaUzation constant and a G (—00,00). 

With constant (3, the Gaussian PDF (20) corresponds to the non-intermittent 
K41 picture of fully developed turbulence, and formally agrees with the experi- 
mental statistics of components of velocity increments in time for the integral time 
scale. However, it fails to describe observed Reynolds-number-depcndent stretched 
exponential tails of the experimental acceleration PDF shown in Fig. 1 which cor- 
respond to anomalously high probabilities for the tracer particle to have extremely 
high accelerations, bursts with dozens of rms acceleration. 

The function /(/?) entering Eq. (16) is a PDF arising from the assumption 
that /? is a random parameter with prescribed external statistics. This is the main 
point of the approach, which extends the usual Gaussian picture. Evidently, the 
characteristic time of variation of the parameter /3 should be sufficiently large to 
justify approximation that the resulting PDF (20) is used with this independent 
randomized parameter. Two well separated time scales in the Lagrangian velocity 
increment autocorrelation have been established both by experiments and DNS.^ 
The large time scale has been found of the order of the Lagrangian integral time 
scale and corresponds to a magnitude part of Lagrangian velocity increments. 

The model (17) belongs to a class of stochastic single-particle models of La- 
grangian turbulence and deals with an evolution of the acceleration in time which 
in accord to the Navier-Stokes equation is driven by time derivative of the rhs of 
Eq. (2). This type of modeling corresponds to the well-known universality^'^^'^^ in 
statistically homogeneous and isotropic turbulence. This is in an agreement with the 
observed temporally irregular character of the Lagrangian velocity and acceleration 
of a tracer particle in high-Re turbulent flows. 

In a physical context, an essential fluid-particle dynamics in the developed tur- 
bulent flow is described here in terms of a generalized Brownian-likc motion, a 
stochastic particle approach, taking the particle acceleration (3) as the dynamical 
variable. Such models are generally based upon a hierarchy of characteristic time 
scales in the system and naturally employ a one-point statistical description using 
Langevin-type equation (a stochastic differential equation of first order) for the dy- 
namical variable, or the associated Fokker-Planck equation (a partial differential 
equation) for one-point probability density function of the variable. 

With the choice of (5-corrclatcd noises such Langevin-type models fall into the 
class of Markovian models (no memory effects) allowing well established Fokker- 
Planck approximation. Consideration of flnite-time correlated noises and the asso- 
ciated memory effects requires a deeper analysis which should be made separately in 
each particular case. The evolution equations are formulated and solved in the La- 
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grangian framework, in a purely temporal treatment, with fluctuations being treated 
along a particle trajectory. Fokker-Planck equation with memory term for joint La- 
grangian single-particle PDF of velocity increments has been studied recently by 

Friedrich.^"' 

Approximation of a short-time correlated noise by the zero-time correlated one 
is usually made due to the timescale hierarchy emerging from the general physical 
analysis of the system and experimental data. Under the stationarity condition, 
one can try to solve the Fokker-Planck equation to find stationary PDF of the 
acceleration, P{a). This function as well as the associated statistical moments can 
then be compared with the experimental data on acceleration statistics. The Fokker- 
Planck approximation allows one to make a link between the dynamics and the 
statistical approach. In the case when stationary probability distribution can be 
found exactly one can make a further analysis without a dynamical reference, yet 
having a possibility to extract stationary time correlators. 

In contrast to the usual Brownian-like motion, the fluid-particle acceleration 
does not merely follow a random walk with a complete self-similarity at all scales. 
It was found to reveal a different, multiscale self-similarity, which can be seen from 
wide tails of a non-Gaussian distribution of the experimental data shown in Fig. 1 . 
This requires a consideration of some specific Langevin-type equations, which may 
include nonlinear terms, e.g. to account for turbulent viscosity effects, and an ex- 
tension of the usual properties of model forces, additive and multiplicative noises. 

Specifically, the class of models represented by Eqs. (16)-(20) is featured by 
consideration of the acceleration evolution driven by the "forces" characterized by 
fluctuating drift coefficient 7 (or fluctuating intensity of multiplicative noise in a 
more general case) and/or fluctuating intensity cr^ of the additive noise. This was 
found to imply stationary distributions of the acceleration (or velocity increments 
in time, for finite time lags) of a non-Gaussian form with wide tails which are a 
classical signature of the turbulence intermittency. Earlier work on such type of 
models arc due to Castaing. Gagne, and Hopfingcr,^"' referred to as the Castaing 
model, in wliicli a log-normal distribution of fluctuating variance of intermittent 
variable was used without reference to a stochastic dynamical equation. 

The difference from the well known class of stochastic models with Gaussian- 
white additive and multiplicative noises which are also known to imply stationary 
distributions with wide tails is that one supposes that intensities of the noises are 
not constant but fluctuate at a large time scale. We refer to the models with such 
random intensities of noises as RIN models. 

This class of models introduces a two-time-scale dynamics, one associated with 5- 
correlation of noises (modeling the smallest time scale under consideration) and the 
other associated with variations of intensities of the noises, their possible coupling 
to each other, and other parameters assumed to occur at large time scales, up to 
the Lagrangian integral time scale. Prom a general point of view, one can assume 
a hierarchy of a number of characteristic time scales. ^'^^ However, as a flrst step 
one simplifies the consideration in order to make it more analytically tractable, in 
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accord to the presence of two characteristic time scales in the Kolmogorov 1941 
picture of fully developed turbulence. 

In the approximation of two time-scales, one can start with a Langcvin-type 
equation, derive the associated Fokker-Planck equation in Stratonovich or Ito for- 
mulations, and try to find a stationary solution of the Fokker-Planck equation, in 
which slowly fluctuating parameters are taken to be fixed. As the next step, one eval- 
uates stochastic expectation of the resulting conditional PDF over the parameters 
with some distributions assigned to them. By this way one can obtain a stationary 
marginal PDF as the main prediction of the model. 

The dynamical model (17) represents a particular simple one-dimensional RIN 
model characterized by the presence of additive noise (a short time scale) and the 
fluctuating composite parameter /? = (a long time scale), where 7 is simply 

kinetic coefficient (a multiplicative noise is not present explicitly) and ct^ is the 
additive noise intensity. This model is, of course, far from being a full model of 
the essential dynamics of fluid particle in the developed turbulence regime. It is 
a theoretical challenge to make a link between the Navier-Stokes equation and 
surrogate one-dimensional Langevin models for acceleration such as that defined by 
Eq. (17). 

The averaging (16) of the Gaussian distribution (20) over randomly distributed 
positive (3, an evaluation of the stochastic expectation, was found to be a simple 
ad hoc procedure to obtain observable predictions, with one free parameter, which 
meets experimental statistical data on the acceleration of tracer particle. One can 
think of this as the averaging over a large time span for one tracer particle, or as 
the averaging over an ensemble of tracer particles, moving in the three-dimensional 
flow characterized by domains with different values of /3 randomly distributed in 
space. 

In a physical context one would like to know processes underlying the random 

character of the model parameter /3. Due to the definition the random character of 
/3 is attributed to a random character of the drift coefficient 7 and/or the additive 
noise intensity a^. In contrast to the usual Brownian motion, in which medium is 
treated thermodynamically in an equilibrium state and therefore parameters en- 
tering dynamical equation are taken constant, the system under consideration is 
characterized by extreme fluctuations and presence of coherent structures that nat- 
urally suggest fluctuating character of some of the model parameters. 

The distribution of (3 is not fixed uniquely by the theory so that a judicious choice 
of f{(3) makes a problem in the RIN model (16)-(20). Below we briefly consider 
three specific models characterized by different prescriptions for distribution of the 
parameter /?, and compare them to the experimental data. 
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3.1.1. The underlying x-square distribution 

With the r (x-square) distribution of /3 of order n (n = 1,2,3,.. .), 



1 / ^ \ 



Wo 



(21) 



the resulting marginal probability density function (16) with P{a\(3) given by the 
Gaussian (20) is found in the form (cf. Ref. 33) 

C 

^^"^ = (a2+n//3o)("+i)/2' ^^^^ 

where C is normalization constant. With n = 3 (/3o = 3 for a unit variance) one 
obtains the normalized marginal distribution in the following simple form: 

This is the prediction of the x-squarc model with the Tsallis entropic index taken 
to be q = (n + 3)/(n + 1) = 3/2 due to the theoretical argument that the num- 
ber of independent random variables at Kolmogorov scale is n = 3 for the three- 
dimensional flow.^^ One can see that the resulting marginal distribution is char- 
acterized by power-law tails that a priori lead to divergent fourth-order moment 
{a^) = JZa^P{a)da. 

A Gaussian truncation of the power-law tails naturally arises under the assump- 
tion that the parameter (3 contains a non-fluctuating part, which can be separated 
out as follows: fi/2 ^ -|- (3/2.'^^ Here, ac is a free parameter measuring the 
non-fluctuating part. This leads to essentially modified marginal distribution 

Cexp[-aVag] 
(a2+n//3o)("+i)/2' ^ ' 

where C is normalization constant and ttc can be used for a fitting. This distribu- 
tion at the fitted value = 39.0 is in a good agreement with the experimental 
probability density function. ^'^ 

Note that at — > oo (no constant part) the model (24) covers the model 
(22). Within the framework of Tsallis noncxtcnsivc statistics, the parameter q — 1 
measures a variance of fiuctuations. For g ^ 1 (no fiuctuations) , Eq. (24) reduces to 
a Gaussian distribution, which meets the experimental data for temporal velocity 
increments at the integral time scale. 

A comparison of the x-square model (23) and x-^quare Gaussian model (24) 
with the experimental data is shown in Figs. 4 and 5. One can see that both the 
distributions follow the experimental P{a) to a good accuracy, although the tails 
of the x-square model distribution depart from the experimental curve at large \a\. 
A major difference is seen from the contribution to fourth-order moment, a^P(a), 
shown in Fig. 5. The x-square model yields a qualitatively unsatisfactory behavior 
indicating a divergency of the predicted fourth-order moment (a*). In contrast, 
the x-square Gaussian model is in a good qualitative agreement with the data. 
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Fig. 4. Lagrangian acceleration probability density function P{a). Dots: experimental data for 
the Rx = 690 flow by Crawford, Mordant, Bodenschatz, and Reynolds.^ Dashed line: the stretched 
exponential fit (4). Dot-dashed line: Beck x-square model (23), q = 3/2. Solid line: the x-square 
Gaussian model (24), ac = 39.0, C = 0.637. a denotes acceleration normalized to unit variance. 

reproducing them well at small and large acceleration values although quantitatively 
it deviates at intermediate acceleration values and gives the flatness value F ~ 46.1 
for flc = 39.0, as compared to the flatness value given by Eq. (5). 

3.1.2. The underlying log-normal distribution 
With the log-normal distribution of f3, 
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Fig. 5. The contribution to fourth-order moment, a*P(a). Top panel: a hnear plot, bottom panel: 
a log-log plot. Same notation as in Fig. 4. 



the resulting marginal PDF (16) with P(a|, 
cently proposed to be'^'' 



P(a) = 



1 

27rs 



d/3 (3 



-1/2 



exp 



given by the Gaussian (20) was re- 
(lnA)2" 



2s2 



(26) 



where the only free parameter s can be used for a fitting, or derived from theoretical 
arguments, = 3 (m = exp[s^/2] to provide unit variance). This distribution is 
shown in Fig. 6 and was found to be in a good agreement with the Lagrangian 
experimental data by La Porta et al.,^ the new data by Crawford et al.,^ Mordant 
et al.,'^ and DNS of the Navier-Stokes equation by Kraichnan and Gotoh.^ 

However, the central part of the distribution shown in Fig. 6 reveals much greater 
inaccuracy of the log-normal model (P(0) ~ 1.23) as compared with that of both the 
X-square and ^-square Gaussian models (-P(O) ~ 0.65) which are almost not distin- 
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Fig. 6. Lagrangian acceleration probability density function P{a). Dots: experimental data for the 
Rx = 690 flow by Crawford, Mordant, Bodenschatz, and Reynolds.^ Dashed line: the stretched 
exponential fit (4). Dot-dashed line: Beck log-normal model (26), s = 3.0. Solid line: Castaing 
log-normal model (27), so = 0.625. a denotes acceleration normalized to unit variance. 



guishable in the region |a|/(a^)^/^ < 4 (Fig. 4); see also recent work by Gotoh and 
Kraichnan.^ This is the main faihng of the log-normal model (26) for = 3.0 al- 
though the predicted distribution follows the measured low-probability tails, which 
are related to acceleration bursts, to a good accuracy. The core of the experimen- 
tal curve (4) (P(0) ~ 0.73) contains most weight of the experimental distribution 
and is the most accurate part of it, with the relative uncertainty of about 3% for 
|a|/(a2)i/2 < 10 and more than 40% for \a\/{a^)^/^ > 40.^ 

The distribution (26) is characterized by a bit bigger flatness value, F = 
3exp[s^] ~ 60.3 for s'^ = 3, as compared to the flatness value (5) which is neverthe- 
less acceptable from the experimental point of view. The peaks of the contribution 
to fourth-order moment shown in Fig. 7 do not match that of the experimental 
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Fig. 7. The contribution to fourth-order moment, a*P(a). Top panel: a hnear plot, bottom panel: 
a log-log plot. Same notation as in Fig. 6. 

curve for the R\ — 690 flow. 

One naturally expects that a better correspondence to the experiment could 
be achieved by an accounting for small scale interactions via turbulent viscosity 
[certain nonlinearity in the first term of the rhs of Eq. (17)] as it implies a damping 
of large events, i.e. less pronounced enhancement of the tails of P{a). 

It should be noted that the idea to describe turbulence intermittency via averag- 
ing of the Gaussian distribution over lognormally distributed variance of some inter- 
mittent variable was proposed long time ago by Castaing, Gagne, and Hopfinger,^"^ 

1 

p{x) = — / dee-'^exp 

where a; is a variable under study. Below, we apply this model to the Lagrangian 
acceleration, x = a. 



2sl 



(27) 
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In technical terms, the difference from the Castaing log-normal model is that 
in Eq. (26) the inverse square of the variance, (3 = 9~'^, is taken to be lognormally 
distributed. In essence, the models (26) and (27) are of the same type, with differ- 
ent parameters assumed to be fluctuating at large time scale and hence different 
resulting marginal distributions. 

One can check that the change of variable, 6 = in Eq. (27) leads to the 

density function different from that given by Eq. (26), 

/•oo 

P{x) ~ / ci/3/3^3/2gxp 
Jo 

where we have denoted mi = thq ^ and si = 2so. Therefore, the distributions 
(26) and (27) are indeed not equivalent to each other, both being of a stretched 
exponential form. 

As to a comparison of the fits, we found that the fit of the Castaing log-normal 
model (27) for the acceleration, with the fitted value .so = 0.625 (too = exp[sQ/2] for 
unit variance), is of a considerably lesser quality as one can see from Fig. 6 and, more 
clearly, from Fig. 7. Positions of the peaks of a^P{a) are approximately the same 
for both the models, namely, |a|/(a^)^/^ ~ 8.0 as compared to |a|/(a^)^/^ ~ 10.2 
for the experimental curve. 

Notice that the existence and positions of the peaks of a'^P{a) reflect some 
characteristic property of the Lagrangian particle dynamics. The value |a|/(a^)^/^ ~ 
10.2 possibly separates different mechanisms underlying stochastic motion of the 
particle. One therefore may expect multiple autocorrelation time scales for a{t). 

We conclude this subsection with the following remark. The Langevin model of 
the type (17), Fokker-Planck approximation of the type (19), and the underlying 
log-normal distribution (25) within the Castaing approach were recently used by 
Hnat, Chapman, and Rowlands^^ to describe intermittency and scaling of the solar 
wind bulk plasma parameters. This model will be reviewed in Sec. 3.2 below. 
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3.1.3. The underlying Gaussian distribution of velocity 

The problem of selecting appropriate distribution of the parameter (3 among possible 
ones was recently addressed in Ref. 36. A specific model based on the assumption 
that the Lagrangian velocity u follows normal distribution with zero mean and 
variance cr^ was developed. The result is that a class of underlying distributions of 
(3 can be encoded in the function (3 = (3{u), and the marginal distribution is found 
to be 



P{a)--C{s) d/3 P{a\f3)exp 



2^2 



du 
d^ 



(29) 



where u{(3) is the inverse function. Note that only an absolute value of u contributes 
to this probability distribution. Particularly, the exponential dependence 



(3{u) = exp[±u]. 



(30) 
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features the log-normal distribution of /? so that Eq. (29) leads to Eq. (26) used in 
Ref. 34, while the x-square distribution of order 1 is recovered when one takes 

(3{u) = V?. (31) 

In general, this model is relevant when ii{u) is a monotonic Borel function of the 
stochastic variable u mapping (—00,00) 9 w to [0,oo) 9 and allows one to rule 
out some ad hoc distributions of 3 as well as to make appropriate generalizations 
of both and log-normal distributions of the parameter /?. 

The possible dynamical foundation of the above model is as follows. The sta- 
tionary distribution (29) with /3(u) = exp[±u] can be associated with the Langevin 
equation of the form^^ 

dta = jF{a) + e'^L{t), (32) 

where we denote u = Tu/2, u follows Gaussian distribution with zero mean, and 
we take 7 = const to simplify the consideration. Here, we adopt a viewpoint that 
statistical properties of the acceleration a are associated with velocity statistics 
due to the Heisenberg-Yaglom scaling (6). Notice also that in Sawford model the 

Langevin equation (10) for a includes velocity u and its variance. 

Below, we outline a relationship of the model (32) to some recent approaches in 
studying the intermittency. 

(i) The form of the last term in Eq. (32), in which w can be viewed as a Gaussian 
process w = ijj{t) independent of the white noise L{t), strikingly resembles that 
involved in the recently developed log- infinitely divisible multifractal random walks 
model by Muzy and Bacry,^^'^° a continuous extension of discrete cascades. This 
model will be reviewed in Sec. 6 below. 

(ii) The driving force amplitude of the form e'^'*\ with the ultraslow decaying 
correlation function of w, 

{u;{t)Lo(t + T))^-\l\n[T/T], (33) 

T < T, in the Langevin-type equation has been recently considered by Mordant et 
al.^ The results of this model have been found in a very good agreement with the 
experimentally observed very slow decay of the autocorrelation of Lagrangian ve- 
locity increment magnitudes l^r^i] for each of the two measured component. Also, 
very slow decay was observed for the cross correlation of magnitudes of compo- 
nents. Both the dynamical correlations were found to vanish only for time scale of 
about Ti, while the autocorrelations of full signed 5rUi decay rapidly (the autocor- 
relation functions cross zero at about 0.06Ti) and the cross correlation function is 
approximately zero. The used time scale r = O.OSTt. It was emphasized that the 
parameter Aq enters both the model autocorrelation functions of velocity increments 
and Lagrangian velocity structure functions. 

This may correspond to the slow e/r decay of the Lagrangian acceleration au- 
tocorrelation predicted by K41 theory. We also note that the fitted value of the 
above intermittency parameter Aq (Aq = 0.115 ± 0.01) is very close to = 1/3, 
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with the value = 3 in Eq. (26) interpreted as the number of independent random 
variables in three-dimensional space at the Kolmogorov scale. If this is not due to 
a coincidence, the intermittency parameter Aq approaches simply the inverse of ef- 
fective space dimension, \q = 1/d, d = 3, for high-Rc turbulent flows. The above 
connection to the MRW model and very slow decay of correlations of absolute val- 
ues of acceleration components indicate relevance of the specific representation (32), 
with very slow varying \u\, in the description of intermittency. In fact, due to the 
experiments^^ the Lagrangian velocity autocorrelation function {u{t)u{t + r)) / (u^) 
decays almost exponentially but very slowly, to vanish only for r > 3Tl, where the 
integral time scale Tl = 2.2 x 10~^ s is two orders of magnitude bigger than the 
Kolmogorov time scale r,, = 2.0 x 10~^ s; R\ = 740 and the mean velocity is about 
10% of the rms velocity. 

(iii) Due to the well-known Kolmogorov power-law relationship between the 
mean energy dissipation rate e and the rms velocity m, e ~ u^/L, the representation 
(32) can be thought of as the result of using the relation ln/3 ~ Ine (statistical law), 
with Ine being normally distributed due to K62 theory.^^ Here, e denotes stochastic 
energy dissipation rate per unit mass. Prom this point of view, one can identify 



where <? is a constant. This means that stochastic dynamics of logarithm of the 
energy dissipation rate is independent, and it influences the acceleration dynamics 
specifically through the intensity of driving stochastic force entering Eq. (32). Sta- 
tionary normal distribution of lo can be in turn derived from the Fokker-Planck equa- 
tion associated with the Langevin equation of a linear form, dtUJ = go+giijj+g2L(t), 
where gi are constants. This equation is in an agreement with the recent results of 
Eulerian (hotwire anemometer) study of the interaction between velocity increments 
and normalized energy dissipation rate by Renner, Peinke, and Friedrich.^^ Partic- 
ularly, they found that an exponential dependence of the diffusion coefficient on 
the logarithmic energy dissipation in the Fokker-Planck equation for the velocity 
increments in space is in a very good agreement with the experimental data. Notice 
that this equation does not imply a logarithmic decay of the Lagrangian correla- 
tion function {u}{t)u){t + r)) proposed by Mordant et al^ This may be attributed 
to the well-known difference between the Eulerian (fixed probe) and Lagrangian 
(trajectory) frameworks. 

(iv) For the choice P{u) = exp['u] corresponding to the log- normal distribution 
of (3, using Eq. (32) one can derive the stationary probability density function of 
the form^^ 



u) = .9 Ine, 



(34) 




(35) 



where g(u) is PDF of u. Hence the joint PDF can be written 



P{a, u) = C{u) exp{ln[g{u)] - 



(36) 
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Such a form of the distribution, containing specificaUy the double exponent, resem- 
bles the "universal" distribution of fluctuations (Gumbel function), 

P{x)=coexp[ci{y-e^)], y = C2{x-cs), (37) 

where Ci are constant, recently considered by Chapman, Rowlands, and Watkins^^ 
(sec also references therein) following the work by PortcUi, Holdsworth, and 
Pinton^^. They used an apparently different approach (not related to a Langevin- 
typc equation) based on the multifractal-type energy cascade and or log-normal 
(K62 theory) underlying distribution of the energy dissipation rate at fixed level. 
They pointed out a good agreement of such P{x) with experimental data, where x 
denotes a fluctuating entity observed in a variety of model correlated systems, such 
as turbulence, forest fires, and sandpiles. The result of this approach meets ours and 
we consider it as an alternative way to derive the characteristic probability measure 
of fluctuations; with g{u) taken to be a (respectively, Gaussian) density function, 
one obtains, up to a pre-exponential factor and constants, P{u) ~ exp[— u — exp(M)] 
(respectively, P(h) ^ exp[— tt^ — exp('u)]). Thus we conclude that the "universal" 
distribution (37) can be derived also within the general framework proposed in 
Ref. 36 that reflects a universal character of the underlying x^ distribution pointed 
out by Beck and Cohen.^^ 

Strong correlation between bursts of Eulerian velocity magnitude and bursts 
of energy dissipation rate has been indicated in the recent work by Pearson et 
al7^ which reports results of DNS of statistically stationary and isotropic slightly 
compressible R\ ~ 219 flow. Unlike to results of some pseudo-spectral methods, 
highly fluctuating character of time series of e, L, u, and R\, where L and u are 
characteristic large length and velocity scales respectively, has been encountered. A 
correct procedure for determination of R\ dependence of = eL/u^ is based on 
accounting for time lapse between the bursts and for the averaging over entire time 
series for which stationarity condition is fulfilled. This leads to elimination of the 
scatter in determining the value of from various recent DNS and experimental 
data. The value — 0.5 has been established as the high-Re asymptotic value, 
which is in agreement with K41 assumption on viscosity independence of the mean 
turbulent kinetic energy dissipation rate for stationary isotropic turbulence. 

In summary, we have presented a class of models (29)-(32) using the basic as- 
sumption that the parameter f3 depends on normally distributed velocity fluctua- 
tions. This class has been found to incorporate the previous RIN models in a unified 
way, with the dependence P{u) required to be a (monotonic) Borel function of the 
stochastic variable u. 

Finally, we note that although successful in describing the observed statistics of 
Lagrangian acceleration, with a few simple hypotheses and one fitting parameter, 
the one-dimensional Langevin RIN models (16)-(19) and (29)-(32) show departures 
from the experimental data and suffer from the lack of physical interpretation in 
the context of the 3D Navier-Stokes turbulence. 
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3.2. Hnat- Chapman-Rowlands model 

A nonlinear Langevin and the associated Fokker-Planck equations obtained by a di- 
rect requirement that the probability distribution satisfies some model-independent 
scaling relation have been recently proposed by Hnat, Chapman, and Rowlands,^'' 
to describe the measured time series of the solar wind bulk plasma parameters. It 
should be emphasized that this approach is related to properties of Fokker-Planck 
equation rather than to those of Langevin-type equation. Nevertheless, we find 
this result relevant to fluid turbulence since it is based on a stochastic dynamical 
framework and leads to the stationary distribution with exponentially truncated 
power-law tails, similar to that implied by the RIN models of Sec. 3.1. 

The Hnat-Chapman-Rowlands model is aimed to describe the observed time 
series of the solar wind bulk plasma parameters and is based on the construction 
of Fokker-Planck equation for which the probability density function obeys the 
following one-parametric model-independent rescaling: 

P{x,t) = t-"°Ps{xt-"°). (38) 

Here, x denotes fiuctuating plasma parameter and ao is the scaling index. Par- 
ticularly, the value ao = 1/2 corresponds to a self-similar (Brownian) walk with 
Gaussian probability density functions at all time scales. The fitted value is differ- 
ent, ao = 0.41, and corresponds to a single non-Gaussian distribution Pg{xs), to 
which the observed distributions of some four plasma parameters collapse under the 
scaling Xs = xt~°'° . 

The Langevin equation of this model assumes only additive noise, and in such 
ansatz it was found to be 

dtx = Di{x) + D2{x)r){t), (39) 
where Di and D2 are of the form 

D,{x) = ,/^xi-"o-V2, (40) 
V ^0 

D2{x) = [60(1 - ao" V2) - aoW-"°' , (41) 

ao and 60 are free parameters and 2Do is intensity of the (J-corrdatcd Gaussian- white 
additive noise ri{t), {rj(t)) — 0. By construction, this specific form of the dynamical 
equation ensures that the corresponding Fokker-Planck equation 

dtP{x, t) = a,[aoa;i-"«" V + 6ox2-"o"'a,P] (42) 

has the general solution P{x, t), which exhibits the scaling (38). The fitted values are 
o-o/bo = 2, 60 = 10, and a^^ = 2.44. The rescaled distribution Pa{xs) corresponding 
to Eq. (42) is characterized by power-law tails truncated by stretched exponential. 
It provides a good fit to the tails of the experimental distribution but diverges at 
the origin Xg ^ 0. 

To sum up, we point out that this diffusion model uses the generalized self- 
similarity principle resembling that used in the Eulerian description of the energy 
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cascade in the developed 3D fluid turbulence and appears to be valid only asymp- 
totically for large values of the variable, with the fitted value of exponent parameter 
being about ao = 0.41. 

3.3. Reynolds stochastic models of Lagrangian acceleration 

The second-order Lagrangian stochastic model of Reynolds''^'''^ which extends Saw- 
ford model (10) is prescribed by the following equation for a component of acceler- 
ation: 

da = -{T^^+r;^^-a-^^'^)adt-TZ\-'udt + ^j2al{TZ^+T^^)TZ^T^^d^. (43) 

Here, is an incremental Wiener process with zero mean and variance dt, 
the energy-containing and dissipative time scales of flow Tl = 2a1^/{Coe) and 
"^rj = Cqv^/"^ / {2aQe^/'^) are defined in terms of the stochastic dissipation rate e, 
universal Lagrangian structure constants ao and Co, the kinematic viscosity v, and 
the velocity variance cr^ = {v?) = . 

Model prediction for the Lagrangian acceleration PDF is due to P{a) = 
P(a|e)/(e)rfe, where one assumes independent zero-mean Gaussian distributions 
for velocity and acceleration with variances cr^ and cr^|^ respectively. This model is 
featured by accounting for fluctuations of e and is consistent with the log-normal 
model reviewed in Sec. 3.1.2. 

In order to calculate P{a) one should determine distribution of the turbulent 
energy dissipation rate e. Following Pope and Chcn,^" the evolution of logarithm 
of the normalized dissipation rate x = ^^{^/^) along a Lagrangian trajectory is 
governed by the stochastic equation 

dx = ~ix - {x))T-'dt + ^2alT^'d^', (44) 

where d^' is independent incremental Wiener process with zero mean and the time 
scale = 2a'^/Co{x}- The distribution of x is thus Gaussian, and its variance was 
approximated by tr^ = 0.354 + 0.289 Ini?;^) in accordance with K62 theory and DNS 
data by Yeung and Popc.^'^ Thus, Reynolds-number effects arc incorporated into 
the model, which is applicable to large time scale. It should be emphasized that the 
introduction of fluctuating x means that the model (43) incorporates both additive 
and multiplicative noises. 

The resulting acceleration flatness factor behaves as F = 1.35R^'^'^, which is in 
agreement with the recent pressure gradient DNS data by Kraichnan and Gotoh^ 
and with lower bound on F set by the experiment.^ The obtained PDF P(o) is in 
agreement with the measured acceleration distribution.-^'^ Also, the model agrees 
well with the observed extended self-similarity of the Lagrangian velocity structure 
functions, the exponential shape of the Lagrangian velocity autocorrelation func- 
tions, and the observed ultraslow correlation of the modulus of the acceleration^'^® 
(see Sec. 3.1.3). 
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Extension of this model to account for the observed dependence of the acceler- 
ation variance on velocity*" will be reviewed in Sec. 5 below. 

One of the important dynamical quantities of a fluid particle motion are tra- 
jectory rotations. Recently developed second-order 3D model by Reynolds and 
Veneziani'*'^ shows that non-zero mean trajectory-rotations are associated with spi- 
ralling trajectories, oscillatory Lagrangian velocity autocorrelation functions, sup- 
pressed rates of turbulent dispersion for given turbulent kinetic energy and mean 
dissipation rate, and skew diffusion; see also work by Borgas, Flesch, and Sawford.''^ 
The 3D stochastic model equation includes terms of the form Cijki^jUkdt, and the 
resulting stationary joint PDF 



P{a,u) = {2TTaiai) ^ exp 







exp 





(45) 



is of a Gaussian form with respect to and Uj. Here, = ct2/(T£,t^) is the con- 
ditional acceleration variance. The conditional mean acceleration {ai\u) = Cijk^jUk 
entering the distribution (45) endows trajectories with a preferred sense of rotation. 

They showed that rotations of the Lagrangian velocity vector produced by such 
model coincide closely with the intense rotations measured in the recent seminal 
experiment by Zeff et alJ^ and are described by Tsallis statistics"^^'^"'^^ to a good 
accuracy. Model predictions for the rotational statistics of the North Atlantic Ocean 
are found to be in close agreement with simulation data produced by the Miami 
Isopycnic-Coordinatc Ocean Model. 

In a recent paper Reynolds''^ constructed phenomenological third-order La- 
grangian stochastic model. The model describes evolution of the material derivative 
of acceleration, the so called hyper-acceleration a = da/dt, in analogy with the 
second-order model. The ID stochastic equation for o, is 

dd = -(T^'+T-i + t^')ddt - {T^\-'+T£%' + T-%')adt 

-T£'T-%'udt + ^J2al^^^{Tz'+T^'+t^')d^, (46) 

da = ddt, du = adt, two timescales Tl and are defined as above, and is third 
timescale related to the hypothesis of finite hypcr-accclcration variance. The hyper- 
acceleration is assumed to be autocorrelated exponentially on the timescale <^Tn, 
in a fully developed turbulent flow. The hyper-acceleration variance cr?|^„ = cr? — 
ct^/ct^ is taken conditional on both the acceleration a and velocity u through their 
variances. 

As in the case of the second-order model, this model applies to stationary homo- 
geneous and isotropic turbulence with Gaussian velocity and acceleration statistics. 
The hyper-acceleration statistics is Gaussian. The ID toy model given by Eq. (46) is 
determined uniquely by the well- mixed condition. The corresponding Lagrangian ve- 
locity autocorrelation function is consistent with the inertial range and dissipation 
range forms of Lagrangian velocity structure functions by Kolmogorov similarity 
theory. 
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The model timescales are related to the variances of velocity, acceleration, and 
hyper-acceleration by 

= '^HTlt^ + TlU + Tr,U)-\ ol = aliT-^T-' +T-%' + T-%'). (47) 

For ^3 = it follows particularly that cr^ = a-2/(Tj,T^) = aQe^/'^u~^/'^ in an agree- 
ment with Heisenberg-Yaglom scaling law (non-intermittent K41 picture) given by 
Eq. (6); e = /L and ao is constant. 

In fact, the third-order model goes beyond Kolmogorov phcnomcnological the- 
ory by introducing third characteristic time-scale t^, in addition to the conventional 
time scales Tl (energy-containing) and r,, (dissipative) . For ^ one thus expects 
deviations from the K41 predictions. Particularly, corrections duo to intcrmittency 
are known to imply dependence of the Kolmogorov constant ao on Reynolds num- 
ber. While for high-Re flows ao is found approximately constant (however, weak 
deviations such as ao — -R^'^^ can not be ruled out), for R\ < 500 there is a clear 
deviation from the K41 scaling of the acceleration variance. -'^ 

Agreement between third-order model predictions and DNS data for Lagrangian 
velocity structure function and Lagrangian acceleration autocorrelation function is 
found significantly better than that obtained with the second-order Lagrangian 
stochastic model (43). The effects of third-order processes were found comparable 
in magnitude to the effects of second-order processes. 

This means that the third-order dynamics at the characteristic time scale is 
essential for a better description of homogeneous isotropic turbulence within the 
framework of such a Lagrangian stochastic modeling approach. Physical interpreta- 
tion of such (and higher-order) processes in the context of 3D Navier-Stokes equation 
and turbulence phenomenology is one of the open problems. 

Anisotropy effects in turbulence within the framework of Lagrangian stochastic 
third-order model has been described by Reynolds, Yeo, and Lee."*^ The simplest 
3D stochastic equation for the hyper-acceleration in homogeneous anisotropic tur- 
bulence is given by 

dhi = {aij + Cij}djdt + {bij - Cikakj)ajdt + dkbkjUjdt + ^J-2{aij + Cy )Ajfed^fe,(48) 

with dtti — didt and dui = Uidt. Here, a^ = —CoCTij — 2aoCQ^{e/i'y^'^6ij, 
bij = -ao{e' /vY/'^Tij, Cijixjk + bjiaik) = -aijXjk, nj = [cr"^]ij, and aij = (u^Uj), 
Xij = {O'iO'j)^ and \ij = (didj) are Lagrangian velocity, acceleration, and hyper- 
acceleration covariances respectively; = 1,2,3 (summation over repeated 
indices is assumed). When the inverse time scales Cij of third-order processes tends 
to infinity the model (48) reduces to a second-order Lagrangian stochastic model. 

Model predictions were compared with the results of DNS acquired for a tur- 
bulent channel fiow with low Reynolds number Rx ~ 30. The third-order model 
is shown to account naturally for the anisotropy of acceleration variances in low- 
Re turbulent flows and for their dependency upon the energy-containing scales of 
motion. The experimental values Co = 6.0 and ao = 5.5 were taken and the hyper- 
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acceleration variances were subsequently chosen to guarantee consistency with the 
DNS data for the acceleration variances. 

As to the high- Re limit, it was argued that if ts/r^ tends to a finite value when 
R\ — > oo then hyper-acceleration statistics may account for the observed anisotropy 
of acceleration variances in high- Re flows. ^ Alternatively, if ta/r^ ^0 for Rx ^ oo 
then one ends up with the isotropic scaling prediction. It was argued that ^ 
implies in general dependence of the acceleration on the energy-containing scales of 
motion through the dependence of Lagrangian acceleration variance upon for 
each component separately. 

Non-universality of parameters in conventional Lagrangian stochastic models, 
when one tries to fit them to experimental and DNS data, was shown to be a con- 
sequence of truncation at either first or second order and not an inherent deficiency 
of the approach in general. Anisotropy of acceleration variance implies "anisotropy" 
of qq (different values of qq for different components) when = 0. This is not nec- 
essarily the case for 0. The additional dynamical degree of freedom provided 
by the third-order model, i.e. inclusion of third-order processes to describe homo- 
geneous turbulence, thus allows one to keep universal character of some parameters 
when fitting model predictions to the experiments and DNS. 

3.4. Laval-Dubrulle-Nazarenko model of small-scale turbulence 

The above one-dimensional Langevin toy models of Lagrangian turbulence consid- 
ered in Sees. 3.1, 3.2, and 3.3 all suffer from the lack of physical interpretation, e.g. 
of short-term dynamics, or small- and large-scale contributions, in the context of 
3D Navier-Stokes equation. 

The crucial point is to make a link between Langevin-type equations and the 
3D Navier-Stokes equation. This includes determination of statistical properties of 
stochastic terms and the functional form of deterministic terms, as well as their 
dependence on the parameters entering the Navier-Stokes equation, justified for 
fully developed turbulence. Also, some extension of the stochastic equation may be 
required to account for dependence of the parameters on Lagrangian velocity, in 
the spirit of RIN approach of Sec. 3.1, and in correspondence to the Navier-Stokes 
equation as the pressure gradient term in the Eulerian framework can be expressed 
in terms of the velocity owing to the incompressibility condition. Strong and nonlo- 
cal character of Lagrangian particle coupling as the result of pressure effects makes 
it difficult to derive theoretically turbulence statistics from the 3D Navier-Stokes 
equation. One is thus left with more or less justified modeling approaches. 

The Navier-Stokes equation based approach to describe statistical properties of 
small scale velocity increments, both in the Eulerian and Lagrangian frameworks, 
was developed in much detail by Laval, Dubrulle, and Nazarenko;**'* see also recent 
paper by Laval, Dubrulle, and McWilliams.*^ 

This approach is based on featuring nonlocal interactions between well separated 
large and small scales — elongated triads — and is referred to as the Rapid Distortion 
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Theory (RDT) approach. Decomposition of velocities into large- and small-scale 
parts was made by introducing a certain spatial filter of a cutoff type. Within the 
framework of this approach, 3D Langevin-type model of small-scale turbulence was 
proposed. 

The main assumption of the Laval-DubruUe-Nazarenko (LDN) approach to the 
3D Navier-Stokes turbulence is to introduce and separate large-scale (L) and small- 
scale (Z) parts in the 3D Navier-Stokes equation and using the Gabor transform 
(Fourier transform in windows)"*'*''''' 

Ui{xn,km,t) = J f{e\xj - a;;.|)e''='"(^'"-^-)wi«,i)d3a;^. (49) 

Here, the parameter e is such that 27r/Z/ <C £ <C 1 and f{x) is a function which 
decreases rapidly in infinity. The "window" to which Fourier transform applies is 
centered at the point Xm and has the size between small and large scales I and L 
respectively. 

This allows to consider analytically small-scale turbulence coupled to large-scale 
terms, i.e. account for an inter-scalc coupling. Such nonlocal interactions were ar- 
gued to be important in understanding intermittency in developed turbulent flows. 
The other, large-scale, part of the equation can be treated separately (and, in princi- 
ple, solved numerically given the forcing and bomidary conditions) since the forcing 
is characterized by presumably narrow range of small wave numbers, and the small 
scales make little effect on it. Small-scale interactions are modeled by a turbulent 
viscosity and were shown numerically to make small contribution to the anomalous 
scaling (intermittency) in the decaying turbulence. Nevertheless, these arc impor- 
tant when fitting model distribution to the experimental data. The 3D LDN model 
of small scale turbulence was used to formulate surrogate ID LDN model, which 
was studied both in the Eulerian and Lagrangian frames.^* 

The ID LDN toy model of the Lagrangian turbulence implies a nonlinear 
Langevin-type equation for a component of small-scale velocity increments in time.^^ 
Such a toy model can also be viewed as a passive scalar in a compressible ID flow. 
For sufficiently small time scale r it corresponds to the acceleration a of fluid par- 
ticle, a = 5ru/T, and is written as"*^'"*^ 

dta={i-Vik^)a + G^. (50) 

This equation corresponds to a Lagrangian description in the scale space, along a 
wave-number packet, defined by the Gabor transform. Here, 

V, = ^Jvl + S2aV/e2 (51) 

stands for the turbulent viscosity introduced to describe small-scale interactions, 
Vq is the kinematic viscosity, i? is a free parameter, k is the wave number \dtk = 
—k^, fc(0) = ko, to model the RDT stretching effect in ID case], ^ and a± are 
multiplicative and additive noises associated with the velocity derivative tensor and 
forcing of small scales by large scales (the energy transfer from large to small scales). 
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respectively. We refer to the model (50) as the ID LDN-type model of Lagrangian 
acceleration in turbulence. 

The entities a± and ^ in Eq. (50) are surrogate versions of the 3D entities related 
to large- and small-scale parts Ui and Ui of the velocity field as follows:^'* 

(52) 

a-± = cr 2 ■ (53) 

(Ji = dj [UiUj - UiUj + UjUi - TJJxh), (54) 

where the hat denotes Gabor transform (49) and the ovcrlinc stands for certain 
spatial cutofi^ retaining a large-scale part. Statistical properties of all the components 
of ^ and <r_L were studied numerically for decaying turbulence and reveal rich and 
complex behavior. 

One can see from Eqs. (52)- (54) that the noises arc related to the velocities and 
their derivatives, and the additive noise a± is associated with interaction terms 
between the large- and small-scale dynamics. One therefore expects that this noise 
may exhibit both the short- and long-time autocorrelations. Physically, this would 
correspond to vortical structures dynamics of which is essentially characterized by 
two well separated time scales. 

This gives support to the idea that the intermittency is caused also by some 
nonlocal interactions including the inertial-range flow modes and not merely by 
small scales. We remark that one would also like to know the role of the dissipative 
scale in this integrated picture. 

Noisy character of the entities (52) and (53) may not be seen as a consequence 
of the Navier-Stokes equation, which does not contain external random forces at the 
characteristic time scale. In the RDT approach, ^ and a± are treated as independent 
stochastic processes entering the small-scale dynamics (50) owing to the fact that 
the large-scale dynamics is weakly affected by small scales (which corresponds to a 
direct energy cascade in 3D flow) and thus can be viewed, in the flrst approximation, 
as a given force of a stochastic character with certain ascribed autocorrelation along 
a particle trajectory. 

Eqs. (50) and (52)-(54) could be used to trace back the origin of multiplicative 
and additive noises entering various surrogate Langevin-type models of the devel- 
oped turbulence, and to provide important information on the dynamics underlying 
the intermittency. 

As a first step, in ID case these noises were modeled in the Lagrangian frame 
by coupled Gaussian-white noises^* inspired by the Kraichnan ensemble used for 
turbulent passive scalar and the Kazantsev-Kraichnan model of turbulent dynamo, 

{a±{t)) = 0, {a±{t)a±{t')) = 2ad{t - t'), (55) 



-^{k-U)-U 
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where D, a, and A are free parameters depending on scale via ko. The positive 
parameters D and a measure intensities of the multipUcative and additive noises 
respectively, while A (to be not confused with the Taylor microscale) measures 
correlation between the noises. 

The representation (55) puts an obvious limitation but is partially justified by 
DNS.^^ The averaging in Eq. (55) is made over ensemble realizations. Zero means 
correspond to isotropy of the stochastic forces along a trajectory. Physically, the 
small scales are thus assumed to be essentially distorted in a stochastic way, as 
a combined effect of much larger scales. The model (55) accounts for short-time 
autocorrelations (approximated by zero- time autocorrelations), with the parameters 
D, a, and A treated here as constants along a particle trajectory. 

We stress that the correlation between the noises ^ and a± defined by Eq. (55) is 
not ad hoc assumption but a consequence of the structure of their 3D counterparts 
(52) and (53) as they contain the same large-scale velocity serving as a unifying 
agent between the noises. 

Partial support of the argument that large scales influence much smaller scales 
is due to the recent study by Pearson et alJ^ of the so called "zeroth law" of 
turbulence. Comparing Lumley's forward cascade model prediction and DNS data 
they argue that some amount of energy is passed to all higher wave numbers, not 
totally to the neighboring wave numbers, at least for low R\ ~ O(IO^). Whether 
this holds for higher-Re flows is however still an open question. 

Stationary solution of the Fokker-Planck equation associated with Eq. (50) and 
the noises (55), 



dtP{a, t) = daiiytk^P) + Dda{adaaP) - Xda{adaP) - \dl{aP) + ad^P, (56) 



where C is a normalization constant and six parameters can be used to fit the 
experimental data. This model specifies the one-dimensional LDN model (50), and 
we refer to this model as the LDN model with (5-correlated noises (dLDN model). 

Thus, one makes a closure by treating the combined effect of large scales, for 
which one has a different dynamical LDN equation that could be in principle solved 
numerically,^^ and nonlocal inter-scale coupling as a pair of given external noises. 
The price of the simplification (55) is that one introduces free parameters a, _D, 
and A to the description. Matching large-scale dynamics to boundary conditions 
deserves a separate study. Despite 3D turbulence is known to be more sensitive to 
the large-scale forcing or boundary conditions, as compared to the 2D one, the used 
simplification (55) is relevant for high-Re flows to some extent**'*'''^ and allows one 
to advance in analytical treatment of the problem. 

It should be stressed that the ID LDN toy model (50) as well as its particular 
case, dLDN model, have several limitations related to the LDN separation of small 



is given by' 



,44 




(57) 
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and large scales allowing to study exclusively nonlocal effects associated with the 
linear process of distortions of small scales by a strain produced by large scales; 
the use of model turbulent viscosity; and one-dimensionality. Applicability of the 
model at Lagrangian integral time scale is limited due to the used separation of 
scales. Formulation of the model to include description at the integral time scale in 
a consistent way is of much interest. Our remark is that one of the possible formal 
ways to accoimt for Reynolds-number effects is to assume that intensities of the 
external noises, which represent large scales, depend on Reynolds number. 

Langevin-typc equation containing both the Gaussian-white multiplicative and 
additive noises was studied in detail by Nakao.^^ The associated Fokker-Planck 
equation was also analyzed. The dLDN model (50) extends Nakao's set up by in- 
corporating two new features: (i) the nonlinearity controlled by B in Eq. (51) and 
(ii) the coupling of the noises controlled by A in Eq. (55). 

It is interesting to note that the RDT approach qualitatively resembles the 
model studied by Kuramoto and Nakao,^* a system of spatially distributed chaotic 
elements driven by a field produced by nonlocal coupling, which is spatially long- 
wave and temporally irregular. Such systems, in which the multiplicative noise is 
the local Lyapunov exponent fluctuating randomly due to the chaotic motion of 
the elements, show power-law correlations, intermittency, and structure functions 
similar to that of the developed fluid turbulence. 

Finally, we note that different strategy to obtain single-particle PDF of La- 
grangian velocity increments from the 3D Navicr-Stokes equation, without referring 
to Langevin-type equations, is due to the Lagrangian PDF method^^ outlined in 
Introduction. 

Before turning to a comparison of the LDN and RIN models which will be made 
in Sec. 4 below, it is worthwhile to outline recent results on RDT approach to 
small-scale turbulence. 

Within the framework of RDT approach, new physical measures of intermit- 
tency such as the mean polarization and the spectral flatness have been introduced 
recently by Dubrulle, Laval, Nazarenko, and Zaboronski^® by using the Gabor trans- 
form of separated large- and small-scale velocities and a method previously devel- 
oped for the kinematic dynamo problem. The resulting equation for the small-scale 
Gabor transformed velocity Um = Um{k^x{t),t), for the fluid-particle trajectory de- 
termined by diXi{t) = Ui, has been found in the following form: 

where aij{t) = VjUi is the large-scale strain. The nonlinear advection term entering 
the 3D Navicr-Stokes equation has been neglected. Hence local interactions among 
small scales are not considered. Nevertheless, Eq. (58) is of much interest since one 
can investigate contribution of nonlocal interactions. 

Gaussian white-noise processes for large-scale strain matrix components (rapid 
stochastic distortions of small scales), the form of which ensures statistical isotropy 



Stochastic models of Lagmngian acceleration of fluid particle in developed turbulence 43 



near a fluid-particle path and the incomprcssibility, 

(o-y (T)cTfc((0)) = il{Sik5ji - d~^Sij5ki)S{T), (59) 

have been taken as a model representation; i. j, . . . = 1, 2, . . . , d. This representation 
in Eq. (58) was used to derive time-evolution equation for the generating function 

Z{k,t) = {exp[Xi\um\'' + X2ul, + Xsu*^]), (60) 

where Ai.2,3 are auxiliary parameters, is constant, and li = 3 for three dimensions. 
Due to the isotropy the final equation depends only on the module of wave number, 
k = \k\. The model (58)-(59) is referred to as Stochastic Distortion Theory (SDT) 
model. 

From the equation for Z an evolution equation for the Gabor- velocity correlators 
of even order can be extracted in a standard way. Inviscid regime and dissipative 
regime approximations of the obtained correlators have been presented. 

Fourth-order wave-number space correlators arc shown to carry essential infor- 
mation on the turbulence statistics and intermittency which is not available from 
the two-point spatial correlators widely used in studying turbulence. Particularly, 
the mean polarization carries information on amplitudes and phase difference of 
Gabor modes. All turbulence wavepackets are shown eventually to become plane 
polarized, — (lyf P), i-e. the turbulence tends to be strongly non-Gaussian. 

Also, the spectral flatness, in the inviscid regime, increases as fc^/^ with the increase 
of wave number k that indicates presence of small-scale intermittency. The effects 
of dissipation have been quantified. 

It was also argued that the log-normal character of turbulence statistics (non- 
linearity in the order of correlator appears to be of a square form in the exponent) 
arises because the strain in SDT model is a multiplicative noise for the velocity, 
which becomes nearly one-dimensional, and the time integrated strain tends to be- 
come a Gaussian process. 

We note that the used model statistics of strain (59) is characterized by S- 
correlations and does not contain long-time fiuctuating terms. Such fiuctuations 
could be accounted for by taking, for example, fl to be long-time correlated stochas- 
tic process, which corresponds to a slow-varying random intensity of the noisy strain 
along a particle trajectory. 

In the following Section, we make a comparison of the RIN model (16)-(19) with 
the LDN-type model (50), as well as its particular case, the dLDN model (57).^^ 

4. Comparison of the simple RIN and LDN-type models 
4.1. A qualitative comparison 

A direct comparison of the Langevin equations (17) and (50) of the two models 
suggests the following evident identifications: 

7F(a) = (^ - zvt/e>, aL = a±. (61) 
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Hence the additive noises can be made identical to each other by putting = a. 
Further, in the case of a linear drift force, F{a) = —a, and constant viscosity, 
ft = vo, we can identiiy the remaining parameters, 'y = I'ok'^ — ^, so that we get 

P = j/a^ = {vok^ - i)/a. (62) 

This relation implies that the parameter /? can be viewed as a stochastic variable 
with a nonzero mean due to the stochastic nature of ^ assumed in the LDN model. 
This is in agreement with the simple RIN model, the defining feature of which is just 
that the fluctuating part of (3 follows some statistical distribution. Such a general 
procedure for obtaining Tsallis-type statistics was suggested in Rcf. 60. 

In the dLDN model (50)- (57), both the additive and multiplicative noises are 
taken S-correlated due to Eq. (55). This is in a sharp contrast to the assumption 
that f3 can be taken constant to derive the stationary solution (20) which is the 
foundation of the simple RIN model. More precisely, the solution in the form (20) 
can be obtained as the lowest-order approximation if /3 is slow varying in time 
as compared to a typical time scale associated with the additive noise L{t) (the 
adiabatic approximation). This suggests that the multiplicative noise ^ should be 
taken as a sufficiently slow varying stochastic variable, to meet the ansatz used in 
RIN models. 

The detailed numerical analysis of the noises^'' for the turbulent flow at rela- 
tively low Reynolds numbers, 57 < Rx < 80, shows that the autocorrelation of the 
multiplicative noise ^ decays much slower (by about one order of magnitude) than 
that of the additive noise a±. Hence the typical time scale at which ^ varies is 
considerably bigger than that of (t±. Also, the cross-correlation between the two 
noises was found to be rather weak, i.e. X <^ D and A <C a, by about two orders of 
magnitude in the longitudinal case, and A = in the transverse case. Altogether this 
allows one to introduce the time-scale hierarchy » Ta and to decouple the noises, 
i.e. to put A = 0, which justifies the adiabatic approximation and gives support to 
one-dimensional RIN models. 

The presence of the long-time correlated amplitude e'^^*^ and the short-time 
correlated directional part L{t) of the stochastic driving force in the Langevin- 
type equation considered by Mordant et al^ also supports the above adiabatic 
approximation (two well separated time scales in the single additive stochastic force, 
in the Lagrangian framework). Also, as established by Hill^'^ for locally isotropic 
turbulence, fluid-particle acceleration correlation is governed by two length scales: 
one arises from the pressure gradient, the other from the viscous force. 

As usual, the 5-correlated noise originates from taking the limit of zero correla- 
tion time in a system with the smallest finite correlation time of the noise. 

On the contrary, in the dLDN model one assumes the approximation of com- 
parable time scales, ~ t^, and retains the coupling parameter A (which relates 
small-scale stretching with vorticity in the 3D case, and is responsible for the skew- 
ness generation along the scale). 
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The use of the constant turbulent viscosity ft = makes a good approximation 
in describing intermittency corrections since both the constant and turbulent vis- 
cosities were found to produce corrections which are of the same level as the DNS 
result. In the physical context, this means that the small scale interactions are 
not of much important in the dynamics underlying the intermittency. This justifies 
the use of the approximation of linear forcing F{a) = —a in simple RIN models. We 
note that this is also in an agreement with both the experimental results for the La- 
grangian velocity autocorrelation function by Mordant, Metz, Michel, and Pinton,^^ 
and the recent experimental Eulerian results for the spatial velocity increments by 
Renner, Peinke, and Priedrich.^* 

Alternatively, one can consider a more general RIN model characterized by the 
presence of Gaussian-white additive and multiplicative noises and fluctuating inten- 
sities of both the noises. This will lead to a model similar to the dLDN model (50) 
in which the noise intensities D and a, and the coupling parameter A are assumed 
to fluctuate at large time scale. 

In summary, we found that the one-dimensional RIN model (16)-(19) can be 
viewed as a particular case of the one-dimensional LDN-type model (50) of turbu- 
lence which is based on the RDT approach by Laval, DubruUe, and Nazarenko.'''* 
It should be stressed that while both the toy models assume introduction of some 
external statistics — the correlator of L{t) and the distribution f{/3) in Eq. (16) and 
the correlators of ^ and a± in Eq. (50) — the LDN-type model is characterized by 
a solid foundation and reveals a rich structure as compared to RIN models. 

In the first approximation, i.e. A = 0, i^t = z^o, and ^ t^, the class of RIN 
models is in a quite good qualitative correspondence with the LDN-type model (50) 
and differs from the specific dLDN model (50)- (57) by the only fact that in the latter 
one assumes ~ and introduces a 5-correlated multiplicative noise. Hence the 
different resulting probability density functions for the acceleration of fluid particle 
in the developed turbulent flow, Eqs. (24)-(26) and (57), respectively. 

4.2. A quantitative comparison 

With the above result of qualitative comparison, we arc led to make a more de- 
tailed, quantitative comparison of the dLDN model (50)- (57) and the simple RIN 
model (16)-(19) with the underlying or log-normal distribution of /3, in order 
to determine which approximation, ~ or ^ Tg-, is better when used to 
describe the Lagrangian statistical properties of the developed turbulent flow. We 
take the recent high precision Lagrangian experimental data^'^ on statistics of fluid 
particle acceleration in the developed turbulent flow as a testbed. Actually we fol- 
low the remark made in Ref. 44 that the (^-approximation of ^ is debatable and the 
performance of such a model should be further examined in the future. 

In Ref. 44, explicit analytic evaluation of the distribution (57) is given for the 
particular case Ut = vq, while the general case is treated in terms of din P {a) /da 
when fitting to the numerical RDT data. In order to make fits to the experimental 
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probability density function P{a) and to the contribution to the fourth-order mo- 
ment, a^P(a), covering wide range of the normahzed acceleration, —60 < a < 60, 
one needs in an analytic or numerical evaluation of the rhs of Eq. (57). To this end, 
we have calculated exactly the integral appearing in the expression (57). Despite 
the integral may look simple the resulting expression is rather complicated. The 
exact result is presented in Appendix A.^^ 

The a-nd log-normal distribution-based probability density functions (24) and 
(26) are both realizations of the RIN model and contain one fitting parameter, Uc 
and s, respectively. The result of comparison of fitting qualities of these functions, 
with Oc = 39.0 and s = 3.0, is that the probability density function (26) provides a 
better fit to the experimental data^ on low-probability tails and on the contribution 
to the kurtosis, which summarizes peakedness of the distribution. However, since the 
integral in Eq. (26) cannot be evaluated analjd;ically we will use the distribution (24), 
which provides a better fit to the central part of the experimental distribution, when 
dealing with analytic expressions. 

The dLDN probability density function (57) contains six parameters which can 
be used for a fitting: the multiplicative noise intensity D, the additive noise intensity 
a, the coupling A between the multiplicative and additive noises, the turbulent 
viscosity parameter B, the parameter uq, and the wave number parameter k. 

The parameter k. For the fitting, we can put k = 1 without loss of generality 
since it can be absorbed by the following redefinition of the parameters and B: 

vok!^ vq, Bk B. (63) 

The parameter a. The structure of the rhs of Eq. (57) is such that only four 

parameters out of five can be used for the fitting. For example, one can put a = 1 
without loss of generality by using the following redefinitions: 

vo/a I/O, B/a B, D/a D, A/a A. (64) 

Alternatively, one can put D = 1 provided the following redefinitions: 

vq/D I/O, B/D B, a/D a, X/D -> A. (65) 

The parameter A. Due to Eq. (A. 5) we have c = —iy/Da — A^, which is imaginary 
for Da > A^ and real for Da < X^. For c = 0, i.e. Da = A^, the integral (A.3) 
is finite since divergent terms in F(c) and F{—c) defined by Eq. (A. 4) cancel each 
other. Since the parameter A measuring the coupling between the noises is assumed 
to be much smaller than both the noise intensities D and a,'*'* we put Da > X^ 
in our subsequent analysis. Moreover, the parameter A responsible for the skewness 
can be set to zero since we will be interested, as a first step, in approximately 
isotropic and homogeneous turbulent flows, for which the experimental distribution 
P(a) exhibits very small or zero skewness.^ 

Thus, we can use three redefined free parameters vq, B, and D for the fitting, 
by putting k = 1, a = 1, and A = 0. However, we shall keep k and a in an explicit 
way in the formulas below, to provide a general representation. 
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We start by considering two important particular cases of the dLDN probability 
density function (57): the constant viscosity Vt = vq and the dominating turbulent 
viscosity Vt = B\a\/k. 

4.2.1. Constant viscosity 

At A = (symmetric case) and B = 0, i.e. constant viscosity Vt = vq, using Eq. (A.l) 
in Eq. (57) we get (cf. Ref. 44) 



where C is normalization constant. This distribution is of a power-law type and one 
can compare it with the result (24), which contains a Gaussian truncation of similar 
power-law tails. 

We note that with the identifications, 



the distribution (66) reproduces that obtained in the context of generalized statistics 

with the underlying distribution.^'^ Particularly, for g = 3/2 (i.e. n = 3 and 
/?o = 3) used there, it follows that D/a= 1 and uok"^ = 3. These values can be used 
as estimates. 

It is highly remarkable to note that the two different approaches yield station- 
ary distribution of exactly the same power-law form, with certain identification of 
the parameters; namely, the Gaussian-white multiplicative and additive noises with 
constant intensities and a linear drift term imply P{a) of the same form as that 
obtained in the RIN model with distributed /3, the ratio of the drift coefficient 
to the intensity of the Gaussian-white additive noise. It follows that the effect of 

distributed (3 mimics the presence of the multiplicative noise, and vice versa, in 
this particular case. 

The power-law distribution (66) can be used to get a good fit of the Lagrangian 
experimental P{a) data for small accelerations, e.g. with the normalized values rang- 
ing from —10 to 10, but in contrast to the Gaussian truncated one (24) it exhibits 
strong deviations for large acceleration magnitudes, and for (1 -I- v^k^ / D)/2 < 2 
leads to a divergent fourth-order moment, which is known to be finite. ^'^^'^^ 

Introducing the noise intensity ratio parameter 



P(a) = C(W + a)-(i+'^«'=V^)/2 



(66) 



D/a = 2{q - 1), (1 + vok^/D)/2 = l/(g - 1), 



(67) 



(68) 



and denoting 



K = -{l + vok'^/D)/2, 



(69) 



we can rewrite the normalized distribution (66) as follows (cf. Ref. 77) 




(70) 



48 A. K. Aringazin and M. I. Mazhitov 



where 2F1 is the hypergeometric function. In accord to the analysis made by 
Nakao/'' for small additive noise intensity, i.e. at b ^ 1, this distribution exhibits 
a pronounced plateau near the origin, and the nth order moments, truncated by 
reflective walls at some fixed |a|, behave as a power of b, 

(a") - b-""''''/^ for n > v^k'^/D, (71) 
(a") ~ 6"" for n < v^k'^/D, (72) 

where n > 0. Thus, the truncated moments behave as 

(a") ~ Go + (73) 

where Go.i are some constants and the function H{n) is zero at v^k"^ = and 
monotonically increases to saturate at n, for large values of i^o^^- It should be 
stressed that such a behavior of the moments for small additive noise intensity is 
not specific to the distribution (70) since it gives divergent moments but arises 
after some truncation of it, for example, by means of reflective walls or nonlinearity. 
Particularly, a truncation of the power-law tails of the distribution naturally arises 
when accounting for the turbulent viscosity to which we turn below. 



4.2.2. Dominating turbulent viscosity 

At A = (symmetric case), for the case of dominating turbulent viscosity, ft = 
B\a\/k, using Eq. (A. 2) we get for positive (upper sign) and negative (lower sign) 
values of a: 

, , Cexp{TBfca/Z)±B/cai/2i:)-3/2arctan[(£)/a)V2al} 
= /p, 2 ^ M/2 ' 

where C is normalization constant. One can see that, as expected, the power-law 
dependence is of a similar form as in Eq. (66) but it is exponentially truncated at 
large \a\ owing to the turbulent viscosity term controlled by the parameter B. This 
distribution is similar to the Gaussian truncated one (24) but the truncation is of 
an exponential type and there is some symmetric enhancement of the tails supplied 
by the arctan term. 

Now we turn to the general symmetric case, which provides a link between the 
two particular cases ft = Uq and vt = B\a\/k considered above. 



4.2.3. The general symmetric case 

At A = (symmetric case), from Eqs. (A.5)-(A.7) we have 

c = —id^, C\ = id\d2, C2 = kd\, (75) 

where we have denoted 

di = yfD{Dkh^f^^B^, d2 = Vd^. (76) 
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Note that c is imaginary and the rhs of Eq. (A. 3) is much simphfied yielding a 
symmetric distribution with respect to a ^ — a. The complex entity C2 defined by 
Eq. (A. 7) may be either real (for Dk^v^ > B'^a) or imaginary (for Dk^v^ < B^a). 
In particular, for the case of constant viscosity, 3> -B^, it is real while for the case 
of dominating turbulent part of the viscosity, B^ ^ Vq, it is imaginary, provided 
that the intensities of noises D and a are of the same order of magnitude. These 
two particular cases lead to different final expressions for the distribution, (66) and 
(74), respectively, obtained above. 

In the general case, using Eqs. (75) and (76) in Eq. (A. 3) after some algebra we 
obtain the following expression for the dLDN probability density function (57), at 
A = 0: 



P{a) 



AD^{B'^Daa'^ + k^{Dkv^ + diVt)"^) 



dlk^{D. 



2 



fedi/(2D^) 



(77) 



where C is normalization constant, di is given in Eq. (76), and ft is defined by 
Eq. (51). 

It can be easily checked that Eq. (77) reduces to Eq. (66) at B = 0, while 

to verify that it reduces to Eq. (74) at Vt = B\a\/k it is required to account for 
that di becomes imaginary, returning back to the logarithmic representation due to 
Eq. (A.4), and the use of the identity (A.8). 

The distribution (77) is characterized by the power-law tails, which arc (i) expo- 
nentially truncated and (ii) enhanced by the power-law part of the numerator, with 
both the efi^ects being solely related to the nonzero turbulent viscosity coefficient B 
responsible for the small-scale dynamics. 

We conclude that to provide an acceptable fit of the dLDN model prediction to 
the Lagrangian experimental data^ small-scale interactions encoded in the turbulent 
viscosity are essential. 

Sample fit of the dLDN probability density function P{a) given by Eq. (77) and 
the corresponding contribution to fourth-order moment are shown in Figs. 8 and 
9, respectively. In the numerical fit, we have put, in accordance to the redefinitions 
(64), the wave number parameter fc = 1 and the additive noise intensity parameter 
a = 1 in Eq. (77) and fitted the remaining three parameters z/q, D, and B. One 
can observe a good agreement with the experimental data. Particularly, the dLDN 
contribution to the kurtosis a^P(a) plotted in Fig. 9 does peak at the same points as 
the experimental curve (positions of the peaks depend mainly on D). The core of the 
dLDN distribution shown in Fig. 8 fits the experimental data to a higher accuracy 
as compared with the log-normal model (26) but yet underestimates that of the 
experimental curve. This departure could be attributed to the approximation of 
5-correlated multiplicative noise used in the dLDN model (see discussion in Sec. 4.1 
above). 

Having the general symmetric form of the dLDN distribution evaluated explic- 
itly, Eq. (77), one can derive the acceleration moments (a"), n = 2,4, . . . The as- 
sociated integrals are not analj^ically tractable and can be evaluated numerically. 
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Fig. 8. Lagrangian acceleration probability density function P{a). Dots: experimental data for 
the Rx = 690 flow by Crawford, Mordant, Bodenschatz, and Reynolds.^ Dashed line: the stretched 
exponential fit (4). Dot-dashed line: Beck log-normal model (26), s = 3.0. Solid line: LDN-type 
model (77), k = 1, a = 1, D = 1.130, B = 0.163, uo = 2.631, C = 1.805. a denotes acceleration 
normalized to unit variance. 

We will consider these below in Sec. 5. 

In the most general case (A 7^ 0) the resulting P{a) is given due to an exponential 
of the exact result (A. 3) which will be considered below in Sec. 5. 

To sum up, we have made an important step forward with the dLDN model 
by having calculated P{a) exactly (see Appendix). We have shown that the dLDN 
model is capable to reproduce the recent Lagrangian experimental data on the 
acceleration statistics to a good accuracy. Particularly, we found that the predicted 
fourth-order moment density function does peak at the same value of acceleration, 
|a|/(a^)^/^ ~ 10.2, as the experimental curve, in contrast to the predictions of the 
other considered stochastic models. The presence of the (5-correlated multiplicative 
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Fig. 9. The contribution to fourth-order moment a'^P{a). Same notation as in Fig. 8. 

noise and the nonlinearity (turbulent viscosity) in the Langevin-type equation was 
found to be of much importance. The considered RIN models provide less but yet 
acceptable accuracy of the low-probability tails although they employ only one free 
parameter, which can be fixed by certain phenomenological arguments, as compared 
to the dLDN model, which contains four free parameters. However, we stress that 
in contrast to the LDN model the considered RIN models have a meager support 
from the turbulence dynamics. 

5. Lagrangian acceleration statistics conditional on the velocity 

In a recent paper Mordant, Crawford, and Bodenschatz"^ reported experimental data 
on the probability density function P{a\u) of the transverse acceleration conditional 
on the same component of Lagrangian velocity. 

Sawford, Yeung, Borgas, Vedula, Porta, Crawford, and Bodenschatz^^ have stud- 
ied acceleration statistics from laboratory measurements and direct numerical simu- 
lations in 3D turbulence at R\ ranging from 38 to 1000. For large the conditional 
acceleration variance was argued to behave as 

{a^\u) - u^ (78) 

obtained to a leading order in the same component u (to be not confused with 
the rms velocity u = [v?)^/"^). This is qualitatively consistent with the stretched 
exponential tails of the unconditional acceleration PDF. The conditional mean rate 
of change of the acceleration derived from the data has been shown consistent with 
the drift term in second-order Lagrangian stochastic models of turbulent transport. 
The correlation between the square of the acceleration and the square of the velocity 
has been shown small but not negligible. 
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In a recent paper Reynolds^" developed a self-consistent second-order stochas- 
tic model with additive noise which accounts for dependence of the Lagrangian 
acceleration covariance matrix Qij = (a,aj|u) on Lagrangian velocity u: 



Here, isotropy of acceleration variances is assumed and the longitudinal and trans- 
verse functions /(luj) and describe covariance of components in acceleration 
parallel and orthogonal to the velocity vector. This model extends the second-order 
Lagrangian stochastic model reviewed in Sec. 3.3 to account for the observed de- 
pendency of Lagrangian acceleration statistics on velocity:^'^'' 



Fitting to the DNS data has been made by the polynomials / = /o+/iu^+/2M^ + 
fau^ and g = go + giu^ + g2U^ + gau^-"^^ The observed dependence of the conditional 
acceleration variance (a^lu) on u was partially understood in terms of Lagrangian 
accelerations induced by vortex tubes within which the vorticity is constant and 
outside which the vorticity vanishes. Dimensional arguments were invoked to derive 
the above third-order polynomial structure of the isotropic covariance matrix as a 
function of squared velocity u^.^'' The parameters /i, /2, gi, and g2 are argued to 
be constants, and /o and go are taken positive. The tails of the resulting P(oi) can 
not be expressed in analytical functions but nevertheless one can show that the tails 
of f (|a|) have a log-normal form. 

The inclusion of such conditional acceleration covariances in the model resulted 
in a significant reduction of the predicted occurrence of small accelerations that 
meets experimental and DNS data for unconditional distributions. The cores of the 
resulting conditional acceleration distributions were found to broaden with increas- 
ing in a qualitative agreement with the experiment, and in general they have 
almost the same shape as the unconditional distribution, in accordance with the 
experimental data.'' 

During the course of derivation of the generalized Fokker-Planck for 
single-particle Lagrangian PDF of velocity increments from the 3D Navier- 
Stokes equation, Friedrich^^ considered the two-point two-time acceleration au- 
tocorrelation conditional on Lagrangian velocity and position, {a{r, t)a{r — 
l,t')\u{xo,t'),x{xo,t'))i=yi^t-tiy Here, Xq is the initial position of a fluid parti- 
cle. This approach and the used closure scheme were outlined in Introduction. The 
conditional dependence was ignored in order to get simple approximation to the dif- 
fusion term. This is consistent with K41 theory. It is of much interest to account for 
such a conditional dependence within the framework of this constitutive approach 
since experimental data and DNS show essential dependence of the Lagrangian 
acceleration variance on the velocity. 



H •11 • 

Qii = [m)-g{\u\)]^+g{\u\)5,i. 



(79) 




(80) 
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Recently, the multifractal approach has been used by Biferale et al."^^ to ob- 
tain acceleration moments conditional on the velocity. Particularly, the multifractal 
prediction 

(a2|u)~M4.57 (81) 

agrees well with the DNS data for large velocity magnitudes. The exponent 4.57 
differs from the value 6 predicted by Sawford et al.'^'^ and is very close to the 
Heisenbcrg-Yaglom scaling exponent vahie 9/2 entering Eq. (6). This indicates that 
the averaging of the conditional acceleration variance (81) over Gaussian distributed 
velocity u is consistent with Heisenberg-Yaglom scaling law (see Eq. (68) and remark 
in Ref. 46). 

The experimental data reveal highly non-Gaussian, stretched exponential char- 
acter of P(a|u), very similar to that of P{a), for fixed u ranging from zero up to 
about three rms velocity^ as opposed to the theoretical result that P{a\u) is a Gaus- 
sian in a due to the simple RIN model (20) for arbitrary (3 = f3{u), or due to the 
more general RIN model (35). Similarity between the experimental P{a\u) and P{a) 
suggests that they share the process underlying the fluctuations. 

Below, wc consider this important problem within the framework of the general 
RIN approach.46.48 

The idea is that stretched exponential form of the tails of observed conditional 
distribution P{a\u) could be assigned solely to small time scales, while the marginal 
probability distribution P{a) is developed from P{a\u) at large time scales, in accord 
to the two-time-scale dynamics. 

This requires some modification in simple RIN models reviewed in Sec. 3.1. The 
sole use of the Gaussian-white additive noise, with fluctuating intensity depending 
on u, and a linear force F{a) = —a with fluctuating 7 = 7(u), is not capable to 
explain the stretching in the observed P{a\u), as it implies only Gaussian conditional 
probability density function P{a\u), for any fixed u. However, it is known that 
accounting for the multiplicative Gaussian- white noise in the drift term of Langevin- 
type equation implies stretched exponential tails. 

Hence one can simply follow the dLDN ansatz as a constitutive model (see 
Sec. 4) using the assumption that the additive noise intensity a appearing in the 
stationary probability distribution P{a\D, a, B, vq) given by Eq. (77) depends on u. 
Also, we will generalize consideration by taking the parameter A to be nonzero and 
depending on u as well. 

The stationary acceleration PDF stemming from the stochastic model (50)-(55) 
has been calculated exactly^^ and due to Eq. (A. 3) is given by 

^ CeM-^tk^/D + F{c, a) + F{-c, a)] 
^ ' (£>a2-2Aa+a)V2(2Bfca + i/tfc2)2SAfc/z?=' > 

for constant parameters. Here, C is normalization constant and vt = ft (a), F{c,a) 
and c are given by Eqs. (A.4)-(A.7). The distribution (82) is characterized by the 
presence of exponential cutoff, complicated power-law dependence, and terms re- 
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sponsible for a skewness of the distribution (asymmetry with respect to a ^ —a). 
While the skewness generation can not be directly applied to our one-dimensional 
case, we do not set A to zero for some reasons to be explained below. 

One way when comparing the model predictions with experiment is to make 
a direct fit of the obtained PDF (82) to the experimental data on unconditional 
acceleration distribution by assuming all the parameters and wave number to be 
constant (sec Sec. 4). 

Particularly, this implies a reduction of the original ID LDN model since wave 
number is taken to be fixed so that the artificial ID compressibility aimed to model 
RDT stretching effect in ID case is not considered. We note that the Lagrangian 
acceleration is usually associated with the dissipative scale, and in the present pa- 
per we do not study dependence of the parameters on the wave number. Such a 
dependence for velocity increments was analyzed in Ref. 44 with the expected re- 
sult that for larger scales the velocity increment PDF tends to a Gaussian form. 
The Gaussian form is reproduced also when D ^ and B 0, i.e. the process 
becomes purely additive with a linear drift term. 

Without loss of generality one can put, in a numerical study, k = 1 and the 
additive noise intensity a = 1 by rcscaling the multiplicative noise intensity D > 0, 
the turbulent viscosity parameter B > 0, the kinematic viscosity Uq > 0, and the 
cross correlation parameter A. The particular cases B = and fo = at A = 0, and 
the general case at A = were studied in detail in Sec. 4. Nonzero A is responsible 
for an asymmetry of the PDF (82) and in 3D picture corresponds to a correlation 
between stretching and vorticity (the energy cascade). Particularly, in the Eulerian 
framework the third-order moment of spatial velocity increment {{6iu)^) was found 
to be proportional to cross-correlation parameter, in accord to a kind of generalized 
Karman-Howarth relationship.'*'* 

However, the approximation based on constant parameters does not allow one 
to consider both the conditional and unconditional acceleration statistics within the 
same model. 

In the following Section, we will consider an extension of the model (50)- (57) with 
the solution (82) by assuming certain model parameters in the obtained stationary 
acceleration PDF to be dependent on random velocity.'** 

This extension is compatible with the 3D LDN approach as ^ and a± depend 
on velocity and contain large-scale quantities due to their definitions (52) and (53). 
Such a functional dependence and the associated longtime fluctuations have been 
ignored when making the simplification (55) with constant parameters D, a, and A. 
We partially restore them. This is the main point of our subsequent consideration, 
and the functional form of the distribution is thus due to Eq. (82) with certain 
parameters being now treated as functions of stochastic velocity u. Observations 
are that the acceleration variance does depend on the same component of velocity. 

We point out that characteristic time of variation of the parameters should be 
suHiciently large to justify the approximation that the resulting PDF (82) is used 
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Fig. 10. A comparison of the experimental unconditional Lagrangian acceleration PDF (dots) 
by Crawford, Mordant, Bodenschatz, and Reynolds^ and the experimental conditional Lagrangian 
acceleration PDF at Lagrangian velocity u = (line) by Mordant, Crawford, and Bodenschatz.^ 
a denotes acceleration normalized to unit variance. 

with independent randomized parameters, P(a|Parameters). Two well separated 
timescales in the Lagrangian velocity increment autocorrelation have been estab- 
lished both by experiments and DNS.^ The large timescale has been found of the 
order of the Lagrangian integral scale and corresponds to a magnitude part that 
is in accord to our assumption that the intensity of noise along the trajectory is 
longtime fluctuating. 

5.1. LDN-type model of the conditional acceleration statistics 

5.1.1. Conditional and unconditional Lagrangian acceleration distributions 

The experimental conditional and unconditional distributions, which we denote for 
brevity by Pexptia\u) and Pexpt(a) respectively, were found to be approximately of 
the same stretched exponential form and both reveal a strong Lagrangian turbu- 
lence intermittency. In Fig. 10 we compare Poxpt(a|0) and Pcxpt(a).^ This similarity 
indicates that they share the same process underlying the intermittency. 

Accordingly, in our previous studies'*^'^^'^^'*^ reviewed in Sec. 4 we used the 
result of a direct fit of the PDF (82) to Pcxpt(a), which was measured with a 
high precision: 3% relative uncertainty for |a|/(a^)^/^ < 10.^''^ We assumed that 
the parameters a and A entering Eq. (82) depend on the amplitude of Lagrangian 
velocity u, while D, B, and are taken to be fixed at the fitted values (fc = 1). 
Theoretically, only a and A may depend on velocity due to Eqs. (52) and (53), while 
the other parameters not. 

An exponential form of a{u) has been proposed in Ref. 46 and was found to 
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be relevant from both K62 phenomenological and experimental points of view. 
Particularly, such a form leads to the log-normal RIN model when u is indepen- 
dent zero-mean Gaussian distributed,^^ and yields the acceleration PDF whose low- 
probability tails are in a good agreement with experiments.'^'*''*® Also, we used an 
exponential form of X{u) so that the conditional acceleration PDF (82) takes the 
form 



Such a form was found to provide good fits of (i) the conditional probability density 
function P{a\u) to Pexpt(«|'w); (h) the conditional acceleration variance {a^\u); and 
(iii) the conditional mean acceleration^* (q^I") varioTis u that meet the experi- 
mental data."^ A brief report on these results is presented in Ref. 82. 

However, self-consistent consideration of the model assumes different strategy:*^ 
P{a\u) should be fitted to Pexpt(«|'") and the marginal PDF computed due to 



where g{u) is PDF of independent random velocity, should reproduce Pcxpt(fl)- The 
marginal distribution (84) corresponds to a convolution of the stationary accelera- 
tion statistics with (independent) random velocity. 

The task is thus to fit a variety of the experimental data, both on the conditional 
and unconditional statistics of acceleration, with a single set of fit parameters. For 
this purpose we use the following natural steps. 

First we fit P{a\u) = P{a\a{u), A(u)) given by Eq. (82) to Pexpt(o|w) assuming 
that the parameters depend on u in an exponential way, 



Hereafter, we use acceleration a and velocity u normalized to unit variances. The 
fit parameter set is D > 0, > 0, i? > 0, Ao, Ua > 0, and > (ao = 1, fc = 1). 
The relations in Eq. (85) mean that the additive noise intensity and the correlation 
between the noises become higher for increasing amplitude of velocity. 

We fit -P(a|0) to Pcxpt(a|0), that excludes the free parameters Ua and ux from 
consideration, by varying D, vq, and i? at ag = 1 and Aq = —0.005. We notice that 
the available conditional statistics Pexpt(a|M) is low for high velocities, the presented 
acceleration range is small, — 14 < a < 14, so that a rather big uncertainty remains 
when determining fit values of the parameters. Changes in shape of Pc^-pt{a\u) with 
u increasing from it = to u = 3.1 are captured independently by the fit parameters 
Ua and u\. The result is shown in Fig. 11. Good overlapping of each theoretical curve 
with data points at all fixed magnitudes of u has been achieved. 

Second we calculate the conditional mean {a\u) and the conditional variance 
{a^\u) and compare them with the experimental data. This decreases uncertainty 
in fit parameter values. The results are shown in Figs. 12 and 13. Note that the 
predicted conditional mean {a\u) as a function of u appears to be very small and 



P{a\u) = P{a\a{u),\{u)). 



(83) 




(84) 



a{u) = aoexp[\u\/ua\, X{u) = Aq exp[|u|/MA]- 



(85) 
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Fig. 11. Theoretical conditional Lagrangian acceleration PDF P{a\ii) (line) and the experimental 
conditional Lagrangian acceleration PDF (dashed line) at velocities n = 0, 0.45, 0.89, 1.3, 1.8, 2.2, 
2.7, 3.1 (from top to bottom, shifted by repeated factor 0.1 for clarity) by Mordant, Crawford, and 
Bodenschatz.^ Lagrangian acceleration and velocity components a and u are normalized to unit 
variances. 



does not match the experiment when the variance is fitted. We wih discuss this 
large discrepancy in Sec. 5.1.2 below. 

Finally we calculate the marginal distribution Pm(a) defined by Eq. (84) by 
using the conditional PDF P(a|a(u), A(m)) and Gaussian distribution of velocity, 



1 



cxp 



u 

T 



(86) 



at fixed a. In the numerical calculation we use a ranging from —100 to 100 with the 
step 0.1. Then we make an interpolation of the obtained points and fit the resulting 
curve to /'cxpt(a)- A noticeable effect of the integration over u with Gaussian g(u) 
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u 

Fig. 12. Theoretical conditional Lagrangian acceleration mean (a\u) (triangles) and the ex- 
perimental conditional Lagrangian acceleration mean by Mordant, Crawford, and Bodenschatz^ 
(squares) as functions of the Lagrangian velocity u normalized to unit variance. 
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u 

Fig. 13. Theoretical conditional Lagrangian acceleration variance (a'^\u) (triangles) and the ex- 
perimental conditional Lagrangian acceleration variance by Mordant, Crawford, and Bodenschatz^ 
(squares) as functions of the Lagrangian velocity normalized to unit variance. 

is a widening of tails of the distribution that meets the experimental data shown in 
Fig. 10. The used integration range for velocity is finite, —20 < -u < 20. 

The fit of Pai{a) to Poxpt(a) strongly decreases the uncertainty but the most 
strict determination of fit values comes due to a comparison of the theoretical con- 
tribution to fourth-order moment, a^P(a), with the experimental data. The results 
are shown in Figs. 14, 15, and 16. Quality of these sample fits is better than in the 
other recent stochastic models reviewed in Ref. 46. In particular, the core of the 
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Fig. 14. Theoretical marginal Lagrangian acceleration probability density function (84) with 
Gaussian distributed Lagrangian velocity (line), experimental data for the Rx = 690 flow (dots) 
by Crawford, Mordant, Bodenschatz, and Reynolds,^ and the stretched exponential fit (4) (dashed 
line). 
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Fig. 15. Lin-lin plot of the central part of the curves of Fig. 14. Same notation as in Fig. 14. 

unconditional distribution reproduces very well that given by the stretched expo- 
nential (4) as one can see from Fig. 15. However, both curves a bit underestimate 
the height at a 0. This means that there is small overestimation in the range of 
medium accelerations, from a — 10 to 30. 

The value Aq = —0.005 has been obtained by adjusting the theoretical curve 
to slightly different heights of the peaks of the observed a'^P{a) shown in Fig. 16. 
Note that the model does not assume the use of ad hoc skewness of the forcing. 
Nonzero cross correlation parameter A naturally results not only in nonzero mean 
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Fig. 16. The contribution to fourth-order moment, a^P(a). Same notation as in Fig. 14. 



acceleration but also in a skewness of both the theoretical distributions P{a\u) and 
Pra{a)- One should extract acceleration components parallel to the velocity vector 
and transverse to it, to verify whether there is skewness in their distributions. 

This skewness could be associated with the Eulerian downscale skewness gen- 
eration, which despite of being small for homogeneous flows is known to be of 
a fundamental character in the inertial range (Kolmogorov four-fifth law), since 
the Eulerian velocity structure function {{6iu)^) was found to be related to cross- 
correlation between two noises. However, we stress that the observed very small 
skewness of acceleration distribution is attributed to the effect of anisotropy of the 
studied flow. How the large-scale anisotropy affects smallest scales of the flow is an 
interesting problem. Our fit made by using nonzero A is of an illustrative character, 
to verify whether it can explain the observed increase of the conditional mean ac- 
celeration with increasing velocity shown in Fig. 12. This issue requires a separate 
study and will be discussed further in Sec. 5.1.2. 

Anisotropic aspects of turbulent statistical fluctuations by using irreducible rep- 
resentations of SO (3) group, which represents a rotational symmetry of the forceless 
3D Navier-Stokes equation, have been reviewed in a recent paper by Biferale and 
Procaccia.^^ They stressed that the anisotropy associated with anisotropic large- 
scale forcing decays upon going to smaller and smaller scales of a high-Re flow, and 
the conflicting experimental measurements on the decay of anisotropy are explained 
and systematized, in an agreement with the S0(3) decomposition theory. 

The following remarks are in order. Our flnding is that the condition Ua < u\ 
provides a convergence of Pni(a)- Also, u\ should not be small to provide assumed 
condition A ^ a at arbitrary u (the cross correlation is small as compared to both 
noise intensities a and D).^'^''^^ We used these criteria when making the fits. 
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The resulting sample fit values arc given by 

D = 2.1, I/O = 5.0, B = 0.35, 
Ao = -0.005, ua = 3.0, ux = 3.0, (87) 

with ao = 1 and A; = 1. The theoretical curves in Figs. 11-16 are shown for this 
sample set of values, which require a further fine tuning. Such a small value of |A| as 
compared to a or is in agreement with that obtained in the LDN direct numerical 
simulations. The calculated flatness factor F = 49.3 of Pm(a) is in agreement with 
the experimental value (5). 

To summarize, the considered Navier-Stokes equation based ID toy model (82)- 
(85) is capable to fit all the available high-precision experimental data on the con- 
ditional and unconditional Lagrangian acceleration statistics^ '^'"^ with the single set 
of parameters (87) to a good accuracy, with an exception being only the conditional 
mean acceleration, which we will consider below. 

5.1.2. Conditional mean acceleration 

As one can sec from Fig. 12, at the values of fit parameters (87) the predicted con- 
ditional mean acceleration {a\u) qualitatively is in agreement but greatly deviates 
from the experimental data. Namely, it is nonzero due to nonzero A and increases 
with the increase of |u| but remains to be much smaller than the experimental data 
even at high values of The conditional mean acceleration is evidently zero for a 
symmetrical distribution (i.e. at A = 0) and should be zero for statistically homo- 
geneous isotropic turbulence when one studies acceleration components aligned to 
fixed directions. The observed departure from zero is thought to reflect large-scale 
anisotropy of the studied Rx = 690 flow. Remarkably, the DNS of (approximately) 
homogeneous isotropic turbulence^ also reveals slightly nonzero mean. In Refs. 7, 
22 a priori statistical equivalence between all the acceleration components in the 
laboratory frame of reference was used and the subsequent averaging over the three 
components gives zero-mean acceleration to a good accuracy. The latter procedure 
is partially justified by the high-Re effect of anisotropy decay when going to smaller 
scales. 

To reduce the discrepancy, one can try the value ux = 1-0 instead of Ux = 3.0 
to provide faster increase of |A| for higher This implies a good fit to the ex- 
perimental conditional mean acceleration (see, e.g. Fig. 2 in Ref. 82) but we found 
an excess asymmetry of P{a\u) at high \u\, with big departure from observations, 
and divergencies when calculating Pyn{a).^^ The reason of the divergency is that 
X{u) at ux = l-O grows faster than a{u) at = 3.0 so that A becomes comparable 
or bigger than a with increasing m, and when Da the function F{c) de- 

fined by Eq. (A. 4) undergoes unbound growth. Thus we conclude that the observed 
nonzero conditional mean acceleration may be due to a persistent effect of the flow 
anisotropy rather than some intrinsic dynamical (cascade) mechanism associated 
with the developed turbulence. It should be emphasized however that some aspects 
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of anisotropic fluctuations in high- Re flows were found to be universal. 

In general one observes a rather small relative increase of the experimental con- 
ditional mean acceleration for higher |u| that eventually reflects a coupling of the 
acceleration to large scales of the studied flow.^'^^"'^ This coupling could be accounted 
for particularly by introducing a correlation between the Lagrangian acceleration 
and velocity. This possibility is of much interest to explore as it may yield the 
deficient increase of {a\u) but it is beyond the scope of the present formalism, 
which assumes an independent velocity statistics. We note also that in contrast to 
the experimental data on the variance (a^jw) the experimental {a\u) exhibits small 
asymmetry with respect to w — > —u (not shown in Fig. 12).^ 

The multiplicative noise intensity D was taken to be independent on the velocity 
u. The effect of variation of D has been considered in Ref. 46 with the qualitative 
result that it does not provide the specific change in shape of P{a\u) observed in 
experiments as shown in Fig. 11. However, a weak dependence of Z) on u can not 
be ruled out. 

5.1.3. Summary 

In summary, the presented ID LDN-type stochastic toy model with the velocity- 
dependent additive noise intensity and cross correlation parameter is shown to cap- 
ture main features of the observed conditional and unconditional Lagrangian accel- 
eration statistics to a good accuracy except for the discrepancy in the conditional 
mean acceleration which can be attributed to certain coupling of the acceleration 
to large scales of the studied R\ = 690 flow.^ 

One of the technical advantages of this model is that one obtains the conditional 
acceleration distribution in an explicit analytical form, Eq. (82), that enables to 
trace back all the produced effects and perform integration over velocity. 

The main result is of course not only good sample fits of a variety of experimental 
data which are important to test performance of the model but also certain advance 
in understanding of the mechanism of Lagrangian intermittency provided by the 
dynamical Laval-Dubrulle-Nazarenko approach to small-scale turbulence. 

The central point is that the LDN toy model has a deductive support from 
Navier-Stokes turbulence. The obtained exact analytic result for the conditional ac- 
celeration distribution and the use of recent high-precision Lagrangian experimental 
data on conditional and unconditional acceleration statistics allow one to make a 
detailed analysis of the mechanism within the adopted framework. 

Effects of large scales and turbulent viscosity have been found of much im- 
portance in steady-state Lagrangian acceleration statistics. The detailed study of 
conditional acceleration statistics reveals a specific model structure of the external 
large-scale dynamics and nonlocal inter-scale coupling for homogeneous high-Re 
flows. Namely, the additive noise a± associated with the downscale energy transfer 
mechanism encodes the main contribution to velocity dependence of the acceleration 
statistics. 
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In general, a cross correlation between model additive and multiplicative noises, 
associated with a correlation between stretching and vorticity in the 3D case, nat- 
urally provides a skewness of distributions and a nonzero mean. Weakness of this 
correlation measured by the parameter A is a theoretical requirement that meets 
the Lagrangian and Eulerian experiments and DNS of homogeneous isotropic turbu- 
lence. It was shown that the observed conditional mean acceleration is not related 
to non-zero A and may be attributed to flow anisotropy effects. In the Eulerian 
frame, the cross correlation between noises is related to the four-fifth Kolmogorov 
law but the effect of skewness generated by A ^ is negligibly small as the result of 
relatively large intensity of the additive noise, which tends to symmetrize accelera- 
tion distributions. This is a dynamical evidence implied by the model rather than 
a result of a priori local isotropy in the spirit of K41 theory. 

The use of exponential dependence of certain noise parameters on statistically 
independent Gaussian distributed Lagrangian velocity has been found appropriate 
to cover new experimental data on conditional statistics and subsequently to transfer 
from the conditional to unconditional acceleration distribution both exhibiting a 
strong Lagrangian intermittency of the flow. Such a dependence is also compatible 
with the log-normal statistics assumed by K62 theory and obtained in the SDT 
model. Gaussian-white multiplicative noise and a longtime correlated intensity of 
Gaussian-white additive noise were both found to make essential contributions to 
intermittent bursts of acceleration. 

It is worthwhile to mention new aspect of stochastic dynamics emerging from 
fluctuating character of the intensity a of additive noise in the presence of the multi- 
plicative noise. The ratio between intensities of these noises, 6^ = D/a, controls the 
character of dynamics of fluid particle. For 6 <C 1 or 6 ^ 1 the acceleration evolution 
is dominated by the influence of additive or multiplicative noise respectively. Since 
b is longtime correlated stochastic parameter, the dynamics can be of a Brownian- 
likc {h <C 1), burst (6 ^ 1), or mixed {b ~ 1) character. The regimes of dynamics 
alternate randomly due to random a and arc characterized by different chances of 
occurrence of bursts of acceleration. Relative weight of each of the regimes in a 
statistically stationary consideration is controlled by the ascribed distribution of a. 

The additive sector of stochastic forces entering the Langevin-type equation 
(50) for the acceleration is of the form e"(*)L(i), where L{t) is Gaussian-white noise 
and u{t) is an independent stationary process with independent increments associ- 
ated with Lagrangian velocity. The stochastic process related to such representation 
has been recently considered by Muzy and Bacry'*^ in the context of multifractal 
random walk*'^'®'* with uncorrelated increments. We will discuss this interesting 
observation^^ in Sec. 6 below. 

6. Multifractal random walk and Lagrangian velocity 

In our previous studies^^'*^'** we took Lagrangian velocity component u in a sta- 
tionary turbulent flow to be Gaussian distributed with zero mean, i.e. u{t) is un- 
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derstood as an independent stationary Gaussian process. For Gaussian distributed 
u the characteristic parameter is the rms velocity u, which completely determines 
statistical properties of the zero-mean random u. The latter is an agreement with 
approximately Gaussian form of the experimental curve for both of two measured 
components of Lagrangian velocity in developed turbulent flows. ^'^^ 

In this Section, we go beyond the Gaussian modeling of the ID process u{t) 
by adopting the following assumption. Lagrangian velocity u{t) is an independent 
stationary stochastically continuous process with independent increments, i.e. the 
increments u{ti)—u{to), u{t2)—u{ti), ^(^3)— ^(^2), ■ • ■, for tk+i > tk are uncorrelated 
and the joint distribution of u{to + t), u{ti + t), u{t2 + t), . . ., for any tk {k = 
0, 1, 2, . . . , n) and any n does not depend on time t. 

In essence, this assumption restates universality of the Lagrangian velocity 
statistics in the inertial range, in accord to K41 similarity hypotheses formulated in 
the Lagrangian framework, and strictly specifies statistical properties of the stochas- 
tic process u{t). 

In general, stationary stochastic processes with independent increments, the sim- 
plest example of which is the usual Brownian motion, is a wide class characterized 
by strong Markovian property (no memory effects) and the so called infinite divis- 
ibility. The most important feature of such processes for the present consideration 
is that one can get a general analytic form of the characteristic function of the 
increments which allows one to determine all statistical properties of the process 
for a given stochastic measure. This enables one to use a broader class of models 
(as compared to the Gaussian modeling) yet having well known analytical tools 
developed by Kolmogorov,^'' Levy,^* and Khinchin.^^ 

By definition, for stationary infinitely divisible processes the characteristic func- 
tion ip[z) of the distribution of u{s + t) — u{s) can always be represented as 
ip{z) = e^^^\ where g{z) is determined by the celebrated Levy-Khinchin formula; see 
Eq. (94) below. For infinitely divisible processes there always exists some function 
iPn{z) such that the relation 

cfiz) = [cpr^izT (88) 

holds for any positive integer n. Indeed, one can use the representation (p„{z) ~ 
exp[i In '^{z)\ , where we fix arg (/^(O) = to provide uniqueness of the representation. 

We start by considering simple example, when the Lagrangian velocity is viewed 
as a stationary continuous process u{t) with independent increments. In many cases, 
without loss of generality one can put the initial value m(0) = 0. This process 
corresponds to the Brownian motion, i.e. u{t) is Gaussian process with the mean 
{u{t)) = 0. We stress that this motion is not due to a molecular structure of the 
medium but is attributed to turbulent fluctuations. This process is stochastically 
continuous (of course, this does not mean that realizations of the process are con- 
tinuous) and infinitely divisible. Indeed, it can be proven that in this case the incre- 
ments u{s+t)—u{s) are a stationary Gaussian process, i.e. u is normally distributed, 
and the homogeneity implies independence on the value of s. Particularly, for the 
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statistical law Ine ~ u we conclude that the stochastic energy dissipation rate e is 
log-normally distributed. This meets K62 log-normal model, in which log-normal 
distribution of e is postulated. 

Let u{t) be a stationary jump process with independent increments, i.e. 

u{t) = const, for t G [0,ti), {tiM), • • • , {tn,T), (89) 

where Q<t\<- ■ ■<tn<T and n is finite number. The velocity is constant at each of the 

finite-time intervals. Again, this process is stochastically continuous and infinitely 
divisible. It can be proven that the jump process is strongly Markovian and the 
number of jumps, j{t), is itself a stochastically continuous stationary Poisson jump 
process with independent increments. With the use of j{t) the process u{t) can be 
represented as 

i(t) 

u{t) = Y^rik, (90) 
fe=i 

where rjk arc independent identically distributed (i.i.d) random variables, and is 
referred to as the generalized Poisson process. 

In particular, if the increments u{s + t) — u{s) are a stationary Poisson process 
with independent increments (the characteristic function is (p{z) = exp[A(e'^ — 
1)]), we conclude that for Ine ~ it the energy dissipation rate e is log-Poisson 
distributed. This meets She-Leveque^° cascade model and extended self-similarity 
inspired model by DubruUe^^ , which are in good agreement with the experimental 
data on scaling exponents of the Eulerian velocity structure functions. 

In general case, a stationary stochastically continuous process with independent 
increments can always be decomposed into the continuous and the jump process 
parts. Therefore, the log-Poisson distribution of e stemming from the jump part 
of the process can receive "corrections" coming from the continuous part of the 
process (log- normal distribution), and vice versa. Hence, the distribution of e arising 
from the general consideration of the stationary stochastically continuous process 
u{t) with independent increments and the statistical law Ine ~ u can be defined 
as a weighted superposition of log-Poisson and log-normal distributions. The weight 
parameter may depend on the scale. 

Muzy and Bacry^^ has recently studied a class of log-infinitely divisible multi- 
fractal random processes of the following form: 

Xr{t) = f e'^-^''Uw{t'), (91) 
Jo 

where Xr{t) = X{t + T) — X{t) is the increment of stochastic process X{t), dw{t) is 
Gaussian white-in-time noise, i.e. w{t) is Wiener process, and w^(t) = InWr/rit)- 
Here, Wr/rit) is a continuous version of independent identically distributed (i.i.d.) 
random variables of the well-known Mandelbrot's discrete multiplicative cascade, 
Mr{t) = WT/T'{t)Mr'{t), in the sense that the scale r takes continuous values in 
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replace of the discrete t„. This corresponds to a continuous extension of the discrete 
cascade. 

The stochastic process uJr{t) can then be represented as a sum of arbitrary 
number of i.i.d. random variables so that by definition LUr{t) is an infinitely divisible 
random process. Note that this process is defined on the scale-time (r, f)-plane, 
with the upper value r = T. The process Ur{t) has the following remarkable scaling 
property: 

UJXr{Xt)=OJr{t) + nr, (92) 

for A < 1, where Sir is an infinitely divisible random variable. 

The multifractal random measure is defined as the small-scale limit, M{dt) = 
limr^o Mridt), of the stochastic measure Mr{dt) = e'^^'^^^'dt. The process Mr{t) = 
Jg* e'^^^*'^rft' is a jump process for certain choice of properties of the characteristic 
function of Wr(0: namely, when the Levy measure has no Gaussian component. 

The small-scale limit, X{t) — limi-^o Xr{t) = B{M{t)), converges and is referred 
to as the Mandelbrot- Taylor subordinated Brownian process. The process X{t) 
represents an example of the multifractal random walk (MRW). 

The MRW can thus be thought of as the usual Brownian motion defined on 
the "multifractal" time r = M{t)^^ Such a viewpoint follows in general from the 
representation of the process X{t) as the stochastic integral of a random function 
f{t) = e'^^*-' over Wiener process w{t). In the case the function f{t) is taken to be 
deterministic one ends up with the usual continuous stationary Gaussian process 
for X[t), i.e. the Brownian motion B(t) defined on usual time t. 

The theory of one-dimensional infinitely divisible stochastic processes developed 
by Kolmogorov, Levy, and Khinchin was applied to the above setup. The result is 
that the multifractal scaling property of the absolute moments oiX{t) can be found 
exactly with the use of the Levy-Khinchin formula for the characteristic function 
ip{z) of the increments of X{t), 



where 



^(^) = es(^\ (93) 



g{z) =imz+ I (e^"" - 1 ~ iz sin x) , (94) 



v{dx) is the canonical Levy measure and m = {X{t)) is the mean. Namely, for the 
case of stationary multifractal measure M{dt) (and finite second-order moment of 
the process) one has*^ 

2P/2„-pr('£±i^ 

{\Xm = ^;^Ty^^P/2^^^"^ (95) 

where Kp = T~''^p^MP) , the temporal scaling exponent is 

C(p) = f-V'(f), (96) 
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and the real convex cumulant-generating function ipiz) is defined by = 9{~i^)^ 
for which without loss of generality one can put '0(1) = 0. The intermittency pa- 
rameter measures nonlinearity of in P and is defined by 

Note also that in Eq. (95) one gets the prefactor in analytically exact form. 

Mordant et al.^ have recently demonstrated that the above MRW formalism in 
constructing of the additive noise in the form e"*^*^L(t), which models the stochas- 
tic force having uncorrelated direction (i.e. L{t) is Gaussian-white) and longtime 
correlation of its magnitude [i.e. u;(i) is Gaussian and ultraslow autocorrelated due 
to Eq. (33)], leads to the Lagrangian scaling exponent C,{p) of the form 

ap) = (]^ + >?)p-^p\ (98) 

C(2) = 1. This corresponds to the log-normal MRW with the generating function 
il>{p) = rnp + where m and n are parameters. Also, the autocorrelation 

functions of logarithm of the amplitude of infinitesimal Lagrangian velocity incre- 
ments in time are in agreement with the experimental data on tracer particle in 
developed turbulent flow. The fit value of the intermittency parameter was found 
to be = 0.115 ± 0.01. We note that this approach does not include a multiplica- 
tive noise sector, such as that represented by the first term in Eq. (50). While the 
ultraslow autocorrelated Lu{t) mimics the presence of multiplicative noise, this phe- 
nomenological approach may result in essential departures from the experimental 
data since the multiplicative noise, which particularly arises naturally in the LDN 
approach reviewed in Sec;. 3.4, is not present explicitly. 

We note that within the same MRW framework one can consider, alternatively, 
log-Poisson MRW, tpip) = [m — sin(ln(5)]p — 7(1 — 6^), which yields a different 
prediction for the Lagrangian scaling exponent 

Cip)=m'p + j{l-SP), (99) 

where m', 7, and 6 are parameters. For ^ — > 1 and 7(ln^)^ — > one recovers 

the log-normal case. The form of this scaling exponent is similar to those of She- 
Leveque^° and DubruUe^^ cascade models corresponding to log-Poisson distribution 
of the normalized energy dissipation rate at spatial scale I, 

<-(p) = (l-A)| + ^^^. (100) 

where A and /3 are parameters, A = /3 = 2/3 corresponds to She-Leveque model 
and C'^(3) = 1; see also work by L'vov and Procaccia.^^ 

The predicted Lagrangian scaling exponents can be fitted to experimental data 
to a good accuracy. The experimental data by Mordant et al.^^ on the relative 
Lagrangian scaling exponents are (C(2) = 1): 

^(p)/<^(2) = 0.56 ± 0.01, 1, 1.34 ±0.02, 1.56 ±0.06, 1.8 ± 0.2 (101) 
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for p = 1,2, 3, 4, 5 respectively. 

Recent DNS results by Biferale et alJ for R\ = 280 flow yield slightly different 
values: C(p)/C(2) = 1.7±0.05, 2.0±0.05, 2.2±0.07 for p = 4,5,6. They mention 
that the range of time delays over which relative scaling occurs is lOr,, to 70t^, 
and in this range anisotropic contributions induced by the large-scale flow appear 
to influence the scaling properties. In the range from to lOr,, all local slopes 
were found to converge around the value 2. This effect is most likely due to the 
capture by very intense small-scale vortical structures, the relative contribution 
of which to the scaling properties is high in the Lagrangian frame, as compared 
to that in the Eulerian one. The influence of dissipative range r < can not 
be associated with these corrections since it tends to increase the value of the 
local slope rather than to decrease it. It was emphasized that the two-time-scale 
dynamics, r G [t^,10t^] and r G [10r,,,Tj,], is a feature of velocity fluctuations 
along Lagrangian trajectories. How to incorporate such dynamical processes in the 
Lagrangian multifractal description and stochastic modelization of particle diffusion 
was mentioned as one of the most challenging open problems arising from their 
analysis. 

In general, the above situation resembles that of fitting to the experimental Eu- 
lerian scaling exponents for which the K62 log-normal model shows an increasing 
departure for large orders p while the She-Leveque and Dubrulle hierarchic models 
are in a good agreement for all experimentally accessible orders p. This indicates 
that log-Poisson models give qualitatively acceptable description of higher inhomo- 
geneities of fully turbulent flow in the Eulerian framework. 

This suggests that within the present model for the longtime behavior the jump 
process approximation is appropriate for qualitative and quantitative description of 
the random walk of the Lagrangian velocity in turbulent flow as compared with the 
continuous process approximation. This can be understood as the fundamental effect 
of flnite size and finite discrete cascade mechanism of the developed turbulence. 

However, both types of processes may contribute to the dynamics, and our main 
conclusion is that one can construct a model which includes them on equal footing 
by adopting the viewpoint that the process u{t) is infinitely divisible with indepen- 
dent increments. Such a process includes both the continuous (Gaussian) and jump 
(Poisson) parts. 

This approach could be used to address the problem pointed out recently by 
Chevillard et al?^ that the Lagrangian (L) and Eulerian (E) singularity spectra 
D{h) can not be both log-normal (parabolic), namely, 

D^{h) = -h+{l + h)D'^{h/{l + h)). (102) 

Whereas the left-hand side of the measured and DNS singularity spectrum curves for 
Lagrangian velocity increments are fitted, the right-hand side of them is reproduced 
by neither the log-normal nor the log-Poisson statistics based models mapped to 
the Lagrangian domain. 

Here we note that this is expected result since intense and weak increments 
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are driven by different processes related to tfie "burst" and "diffusive" dynamical 
regimes. The burst regime corresponds to capturing by very intense vortical struc- 
tures while the diffusive regime describes motion in their incoherent surround. An 
example illustrating the difference is given by the log-normal modeP^ which repro- 
duces low-probability tails of the experimental acceleration distribution to a very 
good accuracy but greatly deviates (overestimates) at the central part of distribu- 
tion which is characterized by low accelerations. 

Within the framework of the ID LDN model it seems to be reasonable to take 
the multiplicative noise ^ to be Poisson distributed and use Gaussian approximation 
for the additive noise (t±. We expect however robustness of the model due to the 
numerical results on the noisy on-off intcrmittcncy model*"^ for which three quite dif- 
ferent types of the multiplicative noise (Gauss-Markov noise produced by Ornstein- 
Uhlenbeck process, dichotomous noise, and chaotic noise produced by Lorenz model) 
were shown to produce the same power-law dependence of the moments for the 
model with reflective walls. 

The above MRW approach gives an additional support to the representation 
e«(*)i(t) and thus to the choice a = e", where u is the magnitude of Lagrangian 
velocity, considered in Sec. 5. Leaving aside the choice of stochastic process for 
u{t) (Gaussian or Poisson, or the weighted superposition of them), the difference 
between MRW approach by Mordant et al.^ and RIN approach to ID LDN-type 
model"*^^'^* is that, in addition to the above additive noise incorporating two well 
separated time scales, the latter model accounts for (i) local interactions via the 
turbulent viscosity, (ii) the multiplicative noise effects, (iii) the coupling of additive 
and multiplicative noises, and (iv) the velocity-dependence of acceleration statistics. 

Appendix A. Exact integrals 

Exact indefinite integrals, up to a constant term which does not depend on a, used 
in calculating the definite integral entering the probability density function (57) are 
given below. 
At vt = vo, 



, -uok^a - Da + X D + i^qP rr. 2 ox , i 

da——-z = — m\Da —2\a+a\ 

Da^-2Xa + a 2D ^ ^ 

Xuok"^ Da — X / , 1 N 

H : arctan— T (A.l) 

DVDa - A2 y/Da - 

At ft = B\a\/k, for positive (upper sign) and negative (lower sign) a, respec- 
tively, 



/ 



^ zpBkd- - D(i + A _ Bka 

" Da^ -2Xa + a ~ ^~D~ 

D^±2BXk, 2 
^ ln[Da^ - 2Xa + a] 



B(Da-2X^)k Da-X , 
± \ , =J=- arctan^ -, (A.2) 
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In the general case, we have obtained a cumbersome expression 
d(^-F^::^^WrV^ = T^-^ HDa -2\a+a] 



Da? -2\a + a D 2 

\n[2Bka + utk^] + F{c) + F{-c), (A.3) 



D2 

where we have denoted, for brevity. 



2D^ 



_c\C2{c — Da + A) 
X (^^(A^ + cA - Da)a + c{Du^k'^ + csz/t))] , (A.4) 



c=-i\/D(x-\^, vt = ^vl + B'^a'^/k'^, (A.5) 
ci = B'^{A\^ + 4cA2 - 3L>aA - cDa) + D^{c + X)u^k^, (A.6) 



C2 = ^B2(2A2 + 2cA - £>a)fc2 + D'^u^k^. (A.7) 
Some useful formulas used in veriiying the limits S — > and D — > are: 

arctana; = 2 (^^(■'- ~ ~ ^^(^ + ^^))' i^-^) 
Um Iln[l + Z)a2]=a2. (a.9) 
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